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SUMMARY 


[luring  the  first  year  of  the  subject  contract,  five 
high-gain,  long-period  ( HGLP ) seismograph  stations  were 
maintained  and  a sixth  station  was  installed  and  maintained. v 
During  the  second  year,  responsibility  f or  mai ntenance  was 
transferred  to  the  Albuquerque  Sei  sinological  Center.  In 
succeeding  years,  Lainont-Doherty  continued  operating  the 
high-gain  station  at  Ogdensburg,  New  Jersey  (OGD).  Research 
completed  under  this  contract  included  investigation  into 
the  characteristics  of  long-period  earth  noise,  the  detection 
capabilities  of  the  high -gain  stations,  digital  signal  en- 
hancement techniques,  and  new  methods  of  distinguishing 
earthquakes  from  underground  explosions,  as  well  as  several 
studies  in  general  seismology  using  high-gain  data. 

'—A  detailed  study  of  earth  noise  at  all  11  HGLP  sites 
found  that  the  spectra  of  earth  noise  at  each  site  is  char- 
actered by  a minimum  in  the  period  range  25  to  45  seconds. 
The  le  1 of  earth  noise  recorded  by  vertical  seismographs 
is  indt  ndent  of  local  meteorological  conditions,  but  at 
shallow  „ites,  the  earth  noise  levels  recorded  by  horizontal 
seismographs  exhibit  a diurnal  cycle  associated  with  daily 
variations  of  atmospheric  turbulence.  The  horizontal  earth 
noise  associated  with  atmospheric  turbulence  is  attenuated 
by  overburden  to  about  10%  of  the  surface  noise  level  at  a 
depth  of  about  150  m.  \ 

An  evaluation  of  the  detection  capabilities  of  the  HGLP 
stations  showed  that  surface  waves  could  be  observed  for  90/, 
of  the  reported  events  of  m^  4.2  or  larger  which  occurred 
within  30°  of  a station.  At  distance  ranges  of  80  to  100  , 
the  detection  threshold  increases  to  about  mu  4.7.  Roughly 
half  of  the  10%  not  observed  are  masked  by  the  coda  of  earliet 
larger  events.  The  thresholds  are  expected  to  be  lowered 
by  about  0.2  units  if  digital  processing  techniques,  such  as 
time-varying  adaptive  filters,  are  applied.  Polarization  and 
azimuthal  filters  are  capable  of  a si gna i -to-noi se  enhance- 
ment of  more  than  6 db  for  signals  from  small  events.  These 
filters  have  proved  to  be  more  effective  than  passive  tech- 
niques or  matched  filtering. 


C A new  surface  wave  magnitude  scale  based  on  the  ampli- 
tude of  higher-mode  Rayleigh  waves  was  developed  and  applied 
to  earthquakes  and  explosions  in  central  Asia.v  several  anom- 
alous events  in  central  Asia,  which  are  characterized  by  un- 
usually low  Ms : «nb  ratios  and  hence  could  be  suspected  to  be 
nuclear  explosions,  are  reclassified  as  earthquakes  when  the 
new,  higher-mode  scale,  Mg,  is  employed.  The  fofcal  depth 
and  mechanism  of  these  anomalous  events  apparently  causes 
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poor  excitation  of  the  fundamental  mode  Rayleigh  wave  with- 
out significantly  affecting  the  amplitude  of  the  higher 
\ mode.  Elimination  of  these  source  related  factors  from 
\tradi  tional  M s : rn  ^ discriminants  is  an  essential  tool  for 
reliable  discrimination  between  earthquakes  and  nuclear 
underground  tests. 

In  other  investigations,  surface  waves  were  used  to 
study  the  focal  mechanism  of  small  magnitude  events,  the 
pulse  distortion  of  body  waves  travelling  non-minimum 
time  paths  was  computed,  and  the  use  of  seismograph  arrays 
in  determining  earth  structure  was  examined. 
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I.  STATION  MAINTENANCE  AND  INSTALLATION 


A.  MARCH,  1971  TO  APRIL,  1972 

Maintenance  was  carried  out  for  the  five  previously  installed  high- 
gain  seismograph  stations  at  Charters  Towers,  Australia  (CTA),  Fairbanks, 
Alaska  ( FLs  K ) , Eilat,  Israel  (EIL),  Chiang  Mai,  Thailand  (CHG),  and  Ogdens 
burg,  New  Jersey  (OGD),  and  a sixth  high-gain  station  at  Toledo,  Spain 
(TLP)  was  installed  and  maintained. 

Preparation  for  the  maintenance  trips  to  CTA,  FBK,  CHG,  and  EIL  began 
in  May  1971.  The  field  teams  involved  in  these  trips  prepared  all  the 
necessary  test  equipment  for  foreign  shipment  and  made  arrangements  for 
their  trips  with  the  different  participating  institutions.  It  was  neces- 
sary to  make  major  modifications  to  the  inverters  and  battery  chargers 
purchased  from  the  Exide  Company  as  part  of  the  uni nterruptable  power 
units  to  be  installed  at  each  site.  Inadequate  wiring  in  the  battery 
chargers  was  corrected  and  relays  and  timers  were  added  to  the  inverters. 
All  necessary  additional  power  cables  and  electrical  equipment  (power 
switches,  fuses,  relays)  were  purchased  and  prepared  for  shipping. 

One  man  from  each  team  spent  as  much  time  as  necessary  to  learn  how 
to  perform  a complete  acceptance  test  on  the  Astrodata  Digital  system. 

One  complete  set  of  spare  digital  cards  and  a tape  deck  were  taken  to  CTA 
and  FBK  to  repair  the  digital  recorders  at  those  sites.  Prior  to  the 
maintenance  trips,  the  digital  systems  at  CTA  and  FBK  were  inoperative 
while  the  header  data  on  the  recorder  at  EIL  was  incorrect.  By  the  end 
of  the  maintenance  trips  to  these  three  stations,  all  three  digital  re- 
corders were  working  properly. 

It  had  become  apparent  that  most  of  the  failures  of  the  digital  re- 
corders incurred  since  the  initial  installation  of  the  high-gain  stations 
were  a result  of  the  station  voltages  being  either  too  low  or  not  suffi- 
ciently regulated,  and  repeated  failures  of  the  station  power  (especially 
at  CTA  and  CliG).  Thus,  in  addition  to  the  uni n terrup tab  1 e power  units 
which  were  installed  at  all  the  sites,  new  power  lines  and  transformers 
were  installed  at  CTA  and  CHG.  The  un i nter ruptab 1 e power  unit,  described 
in  more  detail  in  the  High-Gain  Long-Period  Seismograph  Station  Instru- 
mentation Report  (Lamont,  1971),  consists  of  an  inverter,  battery  charger 
battery  rack,  and  20  six  volt  batteries.  Standard  tests  of  the  seismo- 
graph systems  at  each  site  were  performed.  These  tests  included  redeter- 
mination  of  seismometer  and  galvanometer  free  periods,  damping  character- 
istics, current  sensitivities,  and  absolute  frequency  responses  of  the 
low  and  high-gain  outputs  of  the  three  component  systems.  Copies  of  the 
absolute  frequency  response  determined  at  that  time  were  forwarded  to 
Texas  Instruments  and  the  Sei smol ogi ca 1 Data  Center  in  Asheville,  North 
Carolina,  for  distribution  to  users  of  the  high-gain  data. 

The  maintenance  teams  at  CTA  and  CHG  also  installed  microbarographs 
and  anemometers,  obtained  at  no  cost  to  this  contract.  Data  from  these 
instruments  are  being  used  in  the  study  of  long-period  earth  noise  at 
those  two  sites. 


Failure  of  the  Cipher  tape  decks  in  the  digital  recorders  at  FBK 
on  22  August  1971  and  TLO  on  15  September  1971  necessitated  special 
maintenance  trips  of  approximately  one  week  duration  to  each  site.  The 
FBK  system  was  repaired  on  30  October  1971  and  the  TLO  system  on  5 Nov- 
ember 1971.  Since  these  respective  times  both  digital  systems  have  been 
operational . 

The  digital  recorder  at  CTA  failed  for  the  second  time  on  29  August 
1971.  At  that  time  we  contacted  a digital  technician  based  in  Townsville, 
Australia,  to  service  the  CTA  system.  As  of  15  March  1972,  after  many 
attempts,  he  was  unable  to  repair  the  CTA  recorder.  Thus  arrangements 
were  made  to  ship  a complete  new  digital  unit  from  Lamont  to  CTA. 

On  26  December  1971,  the  digital  technician  from  Townsville,  Australia 
was  sent  to  Chiang  Mai,  Thailand,  to  work  on  the  digital  recorder  at  that 
site,  CHG,  which  had  been  out  of  operation  since  1 Hay  1971.  He  was 
successful  and  CHG  has  been  digitally  operational  since  1 January  1972. 

Lamont  personnel  had,  during  the  time  period  requested,  sent  instruc- 
tions and  replacement  parts  for  many  tests  to  be  carried  out  on  the 
digital  recorder  by  personnel  responsible  for  the  routine  maintenance 
of  the  seismic  equipment  at  CHG.  All  of  these  attempts  were  unsuccessful. 
In  fact,  it  is  now  evident  that  station  operators  at  most  of  the  sites, 
although  proficient  in  their  maintenance  of  the  seismograph  systems, 
cannot  cope  with  the  complexities  of  the  digital  recorders.  For  this 
reason  we  have  attempted  to  contact  qualified  digital  technicians  in  Spain 
Israel,  and  Thailand  who  can  service  the  digital  recorders  in  those  coun- 
tries. He  feel  that  this  is  the  most  efficient  way  to  maintain  these 
systems . 

During  the  time  period  15  July  to  17  August  1971,  installation  of  the 
sixth  high-gain  station  was  completed  at  Toledo,  Spain.  This  station  is 
equipped  with  a digital  recorder  and  an  uni nterruptabl e power  unit.  The 
high-gain  station  is  located  in  a newly  excavated  tunnel  about  4 km  from 
the  Toledo  Geophysical  Observatory.  Personnel  from  the  Toledo  Observa- 
tory, under  the  direction  of  Dr.  Gonzalo  Payo,  are  responsible  for  the 
daily  maintenance  and  record  changing  of  the  high-gain  station.  Peak 
recording  magnifications  obtained  at  this  station  are:  for  the  vertical 

70k  at  40  sec;  E/W  36k  at  47  sec;  N/S  36k  at  40  sec.  The  three  letter 
code  for  this  site  is  TLO. 


B.  HAY,  1972  TO  APRIL,  1973 

The  six  high-gain  seismograph  stations  previously  installed  by  Lamont 
were  maintained  while  carrying  out  an  orderly  transition  of  maintenance 
responsibilities  for  those  stations  to  the  NOAA  group  at  Albuquerque,  Hew 
Mexico.  During  Hay,  June,  and  July  of  1972,  considerable  time  was  spent 
preparing  for  and  carrying  out  combined  inspection-maintenance  trips  to 
the  VLP  stations  at  Charters  Towers,  Australia  (CTA);  Chiang  Hai  , Thailand 
(CHG);  Kipapa,  Hawaii  (KIP);  and  Hatsushiro,  Japan  (MAT).  As  a result  of 
these  trips,  several  recommendations  were  made  for  improvements  such  as: 
the  addition  of  air-conditioners  in  the  recording  rooms  at  CTA  and  CHG 
for  more  efficient  operation  of  the  digital  recorders,  and  Instrument 
repair  techniques  that  should  result  in  greater  continuity  of  the  analog 
data.  In  addition,  personnel  at  each  site  were  thoroughly  advised  of  the 
transition  of  v.is  In  tenant*  respohi  l&  11  i t irs  f r c n Lamont  to  UUAA* 


5 


Poring  the  last  two  weeks  «• 

that""tte  “ludlJs  the  Instrument  tanks,  were  crated  and 
Shi J!ed  ittJef  Si  truck  or  airplane  to  IJMA  at  Albuguerque . The  final 
recordings  from  FBK  are  dated  25  April  1972. 


Seismograms  and  e,.,nftic  tape,  from  all  ‘h^J^tations^r.^ent 
to  Lamont  up  to  1 September  1972  fc r q V c examined  on  a playback  sys- 

?SV lSSlh:55°JS • ;rsSSi.7}S^  ‘ ATfiIrrlS^.9-roup 

iSi5gbihe°^pcsw?or%uai  i ty  «;;i;;;5thcTKps1K52r5S5,I?ea;2nthto°vELA 

inals  forwarded  to  Texas  lnst™mcnts.  in9Ashevi  1 le,  North  Caro- 

Seismological  Center  before  being  fUm  CJ JP^g-qSerSue  has  been  copying 

seismograms  from  any  of  the  stations  at  Lamont. 


C.  HAY,  1973  TO  APRIL,  1974 


Still  KiSisa-is'ss?;  s^* 

special  air  conditioned  room  was  constructed  ®"  "e  J“J“.  J ,nd  h0„ses 
mine  observatory  This  room  Is  12  'J^^nt'g^he  eniironmental  tern. 

per.ture^humidl ty . and  dust  content  1 n ““  UvZitZZ 

ftttjSSft  &B3JB  Mp?sr!5! 

fill  OGD  were  shipped  to  the  Albuquerque  Sei sra-ol ogi cal  Center  in  New 
Mexi co. 


D.  MAY,  1974  TO  APRIL,  1975 


The  hi gh^gain  station  at^Ogdensburg,^tlew^Jersey^(OGD)^was^maintained 

(ASC)  to  replace  the  unit  and  to  recal  brate  jjj  JJj;rJaJy;}80  checkcd 
assistance  of  Lamont  personnel.  The  g J*Jstetirlc  thc  instrument  pres- 
and adjusted  at  this  time.  This  was  * ui  in  a disturbance 

o“fre«h“nth5.rmSl  environment  .S1The  J ffctlSii  gained’,  ns  trumental  ,y 
noisy  for  ^period  of  months  but  gradually  improved. 
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the  ADC  unit  of  the  digital  records  ight  be  at  fault.  Replacement 
parts  were  installed,  but  in  early  December,  the  tape  deck  inexplicably 
began  to  shut  off  intermittently.  It  finally  failed  totally  and  was  not 
restored  to  operation  until  two  technicians  from  ASC  arrived.  At  that 
time,  the  new  solid-state  amplifier  replacements  for  the  old  PTA  were  also 
installed.  One  amplifier  failed  during  installation  and  was  immediately 
replaced.  During  the  first  large  earthquake  following  installation, 
the  magnitude  7.5  event  on  February  2,  1975  in  the  Aleutians,  the  vertical 
and  north/south  solid-state  amplifiers  were  apparently  oversaturated,  and 
ceased  producing  useful  data.  Replacements  were  sent  from  ASC,  but  some 
high-gain  data  was  unavoidably  Ictt.  irior  to  failure  of  the  tol id-state 
amplifiers,  the  records  were  comparable  to  the  old  PTA  records. 

It  is  regretable  that  component  failures  in  the  digital  system  and  the 
high-gain  amplifiers  led  to  a considerable  loss  of  data  in  the  last  year. 
rtuiTai'Iy  2uis>  of  Tost  time  is  attributable  to  digital  tape  fetordtr  failure. 

Both  the  seismograms  and  digital  data  from  OGD  have  been  shipped  to 
ASC  in  New  Mexico. 


E.  MAY,  1975  TO  APRIL,  1976 

The  high  gain  station  at  Ogdensburg,  New  Jersey  has  been  maintained 
in  constant  operation.  Approximately  90  photographic  seismograms  are 
acquired  each  week  in  addition  to  a magnetic  tape  once  per  three  weeks. 
These  data  are  forwarded  to  the  Albuquerque  Se i smol ogi ca 1 Center  and  a 
copy  of  the  tape  is  stored  in  the  Lamont  Data  File.  The  displacement 
outputs  of  the  Geotech  long-period  seismometers  are  recorded  on  2"  per 
hour  analogue  chart  recorders  as  are  the  outputs  from  four  quartz-tube 
s trail. meters . We  are  in  the  process  of  transferring  the  seven  chart 
records  to  a multichannel  3"  per  tiour  recorder. 

The  routine  operation  of  the  observatory  has  required  two  full  time 
technicians  in  the  past.  From  October  1975  we  have  operated  the  station 
with  one  technician  (M.  Connor)  and  have  discontinued  recording  seismic 
data  from  levels  other  than  the  -560  m level  of  the  mine  in  which  the 
observatory  is  located. 

The  digital  recording  of  the  long-period  data  has  been  interrupted 
on  a number  of  occasions  by  electronic  malfunction  of  the  Astrodata 
recorder.  This  is  now  maintained  by  the  Albuquerque  Sei smol ogi cal  Center 
who  have  arranged  on  site  repairs  when  required. 


II.  RESULTS  OF  DATA  ANALYSIS 


A.  EARTH  NOISE 


1 • Earth  Noise  Characteristics  at  High-Gain  Stations 


Until  about  1971  the  bulk  of  the  information  on  earth  noise  at 
periods  longer  than  10  to  20  sec  has  come  from  seismograph  stations 
located  in  the  contiguous  United  States.  The  11  high-gain,  long- 
period  (HGLP ) seismograph  stations  that  were  recently  installed  at 
several  locations  around  the  world  are  providing  important  new 
information  on  the  temporal  and  spectral  behavior  of  earth  noise  at 
periods  longer  than  approximately  10  sec:  An  important  aspect  of 

the  HGLP  stations  is  that  they  are  located  in  a variety  of  geologic 
and  tectonic  settings  (Figure  1).  In  particular,  the  installation  of 
a HGLP  seismograph  station  at  Kipapa,  Hawaii  (kIP),  has  enabled  us 
for  the  first  time  to  examine  the  characteristics  of  long-period 
earth  noise  at  a mid-ocean  site. 


A major  contribution  of  this  study  is  the  significant  increase 
in  the  number  of  long-period  earth  noise  spectra  that  are  available 
from  areas  other  than  the  contiguous  United  States.  The  information 
presented  in  this  study  forms  a comprehensive  set  of  data  on  long- 
period  earth  noise  which  will  be  useful  for  various  seisinol  ogi  cal 
investigations.  For  example  Linde  and  Sachs  (1971)  pointed  out  the 
importance  of  taking  into  account  the  contribution  of  noise  to  the 
long-period  portion  of  the  spectra  of  seismic  waves  in  any  attempt 
to  distinguish  between  differing  theories  of  earthquake  radiation 
(Archambeau,  1968;  Brune,  1970).  As  a second  example,  the  noise 
spectra  presented  in  this  study  pertain  to  detection  and  discrimina- 
tion between  earthquakes  and  underground  nuclear  explosions  in  that 
they  provide  an  estimate  of  minimum  detectable  signals,  especially 
surface  waves  (Rayleigh  and  Love),  as  recorded  by  a singie  instru- 
ment at  a particular  site.  Knowledge  of  the  spatial  and  temporal 
chara'*f  nr' e rs  lonn-neriod  earth  noise  would  aid  in  the  design 
(element' spacing)  of~ seismic  arrays  that  might  further  improve  the 
detection  capabilities  for  surface  waves. 


Digital  data  samples  were  selected  on  the  following  . 

Initially  the  analog  (photographic  paper  or  70  mm  film  chips)  seis  o- 
grams  we^e  examined,  if  they  were  available,  to  obtain four  to  f i ve 
hours  of  .data  during  which  no  reported  earthquakes  or  transient  dis 
turbances  could  be  observed.  The  data  sample  did  not  have  to  b 
continuous  because  the  program  for  computing  the 

c l ■:  n over  predetermined  time  segments.  Next  the  Monthly  List  or  ne 
li mi naryfDe terminations  of  Epicenters  (National  Ocean  Survey,  now 
l.V.  GMlogic”  Survey ) was  checked  for  occurrence. of  unobserved 
Pflrthnuakes  during  the  selected  time  intervals.  Finally,  the  9 
daJa2e?e  inspected  for  transient  disturbances  by  using  a cathode  ray 
iJbe.  Thus?  it  was  assured  that  the  data  selected  for  each  station 
were  free  of  earthquakes  and  spurious  transients. 
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Fig.  1.  Map  of  the  world  showing  the  location  of  the  11  Jiigh-gain,  long-period  stations 
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figure  2. 


Earth  noise  power  spectra  for  the  vortical  seismographs  at  OGD  and 
Z1.P.  The  power  spectra  in  all  the  remaining  figures,  except  figure 
12a  and  b,  have  beer,  set  to  their  absolute  value  and  corrected  for 
instrumental  response.  The  SDi  confidence  interval  indicated  by  the 
error  bars  is  approximately  the  same  for  all  power  spectra  presented. 
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Earth  r.oiso  power  spectra  for 


the  horizontal  seismographs  at  OG'J  and 


Each  power  spectrum  was  estimated  by  time-averagi nq  the  square 
of  the  Fourier  coefficients  of  up  to  twenty  time  segments  of  the 
data  set.  Each  of  the  time  segments  was  1,024  sec  long,  with  the 
exception  of  the  KON  data  segments  which  were  512  sec  long.  The 
spectra  shown  here  are  corrected  for  instrument  response  and  set  to 
the  absolute  level  of  ground  motion  unless  explicitly  noted. 

Deep  Sites  There  are  four  tlGLP  stations  with  more  than  200  m of 
overburden  above  the  seismometer  vaults.  They  included  OGD,  the 
prototype  HGLP  station  at  543  m,  EIL  at  200  m,  KON  at  340  m,  and 
ZLP  at  300  m.  We  will  consider  OGD  and  ZLP  first;  KON  will  be  given 
special  consideration  later. 

Visual  examination  of  all  the  available  70  mm  film  chips  from 
OGD  and  ZLP  shows  that  the  level  of  vertical  earth  noise  at  these 
two  stations  is  independent  of  all  small-scale  meteorological  con- 
ditions of  local  origin.  In  the  case  of  OGD,  nearly  three  years  of 
data  were  examined.  Power  spectra  typical  of  .background  earth  noise 
at  OGD  and  ZLP  are  shown  in  Figure  2.  These  spectra,  which  are  cor- 
rected for  instrument  response,  show  comparable  shapes  and  absolute 
levels  at  periods  longer  than  about  25  sec  and  show  a broad  minimum 
between  25  and  45  sec.  At  periods  shorter  than  about  20  sec,  the 
OGD  spectrum  for  this  particular  time  interval  has  higher  long-period 
(16  to  20  sec)  microseismic  level.  The  fairly  rapid  increase  (12  to 
14  db/octave)  in  the  spectral  levels  at  periods  longer  than  about 
50  sec  is  partially  controlled  by  instrumental  noise  (Savino  et  al . , 
1972;  Fix,  1973).  By  this  statement  we  are  referring  to  the  fact 
that  the  earth  noise  spectrum  is  not  a constant  level  above  the 
instrumental  noise  of  the  seismograph  system.  The  40  sec  power  of 
earth  noise  at  OGD  is  about  20  db  above  the  instrumental  power  while 
at  about  100  sec  the  separation  is  about  10  db.  If  the  instrumental 
noise  of  the  seismograph  system  fluctuates  with  time,  it  is  possible 
that  the  instrumental  noise  may  begin  to  contaminate  the  earth-noise 
spectrum  at  periods  longer  than  about  50  sec. 

Analysis  of  a similar  large  quantity  of  data  from  horizontal 
seismographs  operating  at  OGD  and  ZLP  indicate  that  the  levels  of 
earth  noise  recorded  by  horizontal  seismographs  are  generally  inde- 
pendent of  local  meteorological  conditions.  In  addition  as  indicated 
in  Figure  3 during  ordinary  weather  conditions,  the  shape,  absolute 
levels,  and  the  occurrence  of  a minimum  are  similar  to  those  shown 
in  power  spectra  for  earth  noise  recorded  by  vertical  seismographs 
(Figure  2').  There  is  an  indication  that  the  absolute  levels  of  the 
horizontal  power  spectra  are  lower  than  the  absolute  level  of  the 
vertical  spectra  during  the  same  time  interval. 

Shallow  Sites  There  are  seven  HGLP  s tati  ons  with  less  than  50  m 
of  overburden  above  the  seismometer  vault.  First  to  be  considered 
is  CTA  (30  m of  overburden)  since  both  analog  and  digital  recordings 
from  a single  on-site  microbarograph  are  available  from  this  site, 
permitting  a quantitative  correlation  between  local  meteorological 
conditions  and  earth  noise. 


Typical  examples  of  power  spectra  of  the  earth  noise  at  CTA 
spanning  a two-year  time  interval  are  shown  in  Figures  4a,  b,  and 
c.  In  Figure  4a  note  the  minimum  in  the  spectra  between  about  30 
and  45  sec  and  that  it  is  stable  in  absolute  level  in  the  passband 
to  within  about  6 db  for  the  two-year  time  interval.  The  6 db  vari- 
ation is  barely  above  the  95%  confidence  limit  of  4 dh  for  the  power 
spectra  calculations.  Also  note  that  for  the  same  time  interval  that 
the  power  in  the  long-period  (16-20  sec)  microseism  band  varied  by 
12  to  14  db.  The  minimum  in  the  spectra  <f  horizontal  earth  noise 
tends  to  be  broader  (25  to  45  sec)  and  to  be  less  pronounced  (Figures 
4b  and  c).  In  fact,  the  minimum  is  nearly  filled  in  by  a higher 
noise  level  on  the  horizontal  spectra  for  September  21,  1973.  But 
note  that  the  vertical  spectrum  is  unchanged. 

Nearly  three  years  of  70  mrn  film  chips  were  examined  to  verify 
the  results  given  by  a limited  number  of  power  spectra  spread  over 
a,  two-year  time  period.  The  film  chip  data  and  the  power  spectra 
show  that  earth  noise  recorded  by  the  vertical,  seismograph  does  not 
experience  any  apparent  diurnal  variations  and  the  40  sec  level  of 
earth  noise  recorded  by  the  vertical  seismograph  is  unaffected  by 
variations  of  the  16-18  sec  microseism  level. 

Charters  Towers.  Earth  Noise  as  Correlated  to  Hi crobarograph 
Noise,,  Examination  of  the  70  mm 'film  chips  from  CTA  has  made  it  very 
obvious  that  the  horizontal  earth  noise  undergoes  a very  definite 
diurnal  cycle.  This  cycle  characteristically  goes  from  03:00  •*> 

15:00  noisy  to  17:00  * 02:00  quiet  (local  times).  The  relative 
terms  "noisy"  and  "quiet"  involve  a noise  decrease  of  about  6 db. 
Figure  5 also  shews  that  the  fluctuations  in  the  40  sec  level  of  the 
horizontal  noise  correlate  very  well  with  variations  in  the  40  sec 
level  of  the  noise  recorded  by  the  microbarograph. 

Cross-spectra  between  earth  noise  recorded  by  the  horizontal 
seismographs  and  the  atmospheric  pressure  variations  recorded  by  the 
single  on-site  ini  crobarograph  were  calculated  to  study  the  correla- 
tion between  these  two  types  of  noise.  The  cross-spectra  strongly 
indicate  that  the  variation  of  local  barometric  pressure  is  one  of 
the  principal  sources  of  the  earth  noise  recorded  by  the  horizontal 
seismographs.  McDonald  e t a 1 . (1971)  and  Priestly  ( 1966)  have  ob- 
served that  mi crobarographs  that  do  not  employ  properly  designed 
pipe  array  filters  or  have  no  filter  at  all  (the  situation  of  the 
microbarograph  at  CTA)  are  principally  sensitive  to  the  atmospheric 
turbulence  in  the  boundary  layer  immediately  above  the  earth's 
surface.  Thus,  the  cross-spectra  indicate,  more  specifically,  that 
atmospheric  turbulence  in  the  boundary  layer  immediately  above  the 
earth's  surface  is  a principal  source  of  earth  noise  recorded  by 
horizontal  seismographs. 

Significant  coherence  at  periods  longer  than  about  30  sec  has  been 
observed  during  daylight  hours  between  earth  noise  recorded  by  the 
horizontal  seismographs  and  atmospheric  pressure  variations  recorded 
by  the  microbarograph  (Figure  6b).  In  addition  in  Figure  6a  note  the 
similar  correlation  between  the  earth  noise  recorded  by  the  two  hori- 
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and  of  microbarographic  power  at  CTA  on  seven  days  in  1972.  The 

<°  sec  levels  of  earth  noise  power  are  plotted  at  their  absolute 

values,  the  .,0  sec  level  of  microbareoraphic  power  is  plotted 
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that  the  level  of  vertical  earth  noise  varies  over  a narrow  band 
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band.  The  levels  of  microbarometric  noise  and  horizontal  earth 
noise  are  reasonably  well  correlated. 


COHERENCE  SQUARED  COHERENCE  SQUARED 


SRS5S5SSWB1 


08 


CTA-DAY 

N/S  VS  PRESSURE 


06b 


04k 


0.2  k 


. H . 


J L..1...1 


!00  5040  30  30  15 

PERIOD,  sec 


95  *_  _ 

« . 

i-tl  I 

10 


(C)  (d> 


1.0 

1.0 

pTA- NIGHT 

CTA- NIGHT 

N/S  VS  E/W 

N/S  VS  PRESSURE 

00 

— 

— 

08 

— 

— 

06 

— 

— 

0.6 

— 

— 

• 

• 

04 

— 

— 

04 

• 

-* 

- 

- 

- 

■ f 

• • 

0.2 

nr — • 

• • . 

95 

• 

02 

• 

- • • • • 

_?5%_T7 

• 

f Vr-V'vT  *..i  ' 

• 

• * • M. 

1 • 1 . . I — * 1 J ■ * 

s‘  • 
al  ■ 

100  5040  30  20  15 

1C 

) ° 

100  5040  30  20  15 

10 

PERIOD,  sec 

PERIOD,  sec 

Figure  6.  Squared  coherence  measurements  between  the  two  horizir.tal  seismographs 
(a)  and  between  the  nerth/south  horizontal  seismograph  and  micro'oaro- 
graph  (b)  at  CTA  with  data  samples  taken  during  daylight  hours.  Co- 
herence squared  measurements  between  the  two  horizontal  seismographs 
(c)  and  between  the  north/south  horizontal  seismograph  and  the  micro- 
barograph  (d)  at  CTA  with  data  samples  taken  during  nighttime  hours. 
Note  that  the  squared  coherence  between  the  daylight  measurements 
are  significant  from  about  30  to  100  sec  but  that  the  nighttime 
measurements  display  significantly  less  coherence. 


zontal  seismographs.  During  nighttime  hours  the  correlation 
between  the  eartli  noise  and  atmospheric  pressure  variations  is  sig- 
nificant only  at  periods  greater  than  about  170  sec  (Figure  6d); 
a similar  observation  is  the  case  for  the  correlation  of  the  earth 
noise  recorded  by  the  horizontal  seismographs  (Figure  6c).  A com- 
parison of  the  power  spectra  of  horizontal  earth  noise  from  both 
nighttime  and  daylight  time  intervals  and  of  atmospheric  pressure 
variations  provides  evidence  that  the  majority  of  the  difference 
between  nighttime  and  daylight  power  spectra  of  horizontal  noise 
may  be  attributed  to  differences  in  atmospheric  pressure.  Table  1 
shows  that  if  the  difference  between  nighttime  and  the  daylight 
mi crobarographi c power  spectrum  is  added  to  the  nighttime  power 
spectrum  of  horizontal  earth  noise,  one  obtains  a very  good  approxi- 
mation to  the  daylight  power  spectrum  of  earth  noise.  This  is  shown 
by  the  near  equality  of  the  spectral  differences,  day-night,  of  the 
horizontal  earth  noise  and  spectral  difference  of  the  microbaro- 
graphic  noise.  Note  that  vertical  spectral  differences  are  not  nearly 
equal  to  pressure  differences. 

Sources  of  Earth  Noise  Recorded  by  Vertical  Seismographs  at  all 
HGLP  Sites  and  by  Horizontal  Seismographs  at  Deep  Sites.  Savino  et  al. 
(1972)  showed  that  for  the  HGLP  instrument at i o n a t TJGDT  the  level  of 
earth  noise  at  40  sec  was  approximately  20  db  above  the  instrumental 
noise  levels.  At  200  sec  this  separation  has  decreased  to  about  3 db, 
so  we  have  limited  this  study  to  periods  less  than  200  sec.  Figures 
7a  and  b are  shown  to  emphasize  the  similarity  of  the  power  spectra 
of  vertical  eartli  noise  for  ten  of  the  HGLP  seismograph  stations. 


Despite  observations  that  the  earth  noise  recorded  by  the  vertical 
seismograph  does  not  change  appreciably  with  depth  and  that  the  long- 
period  end  of  power  spectra  of  the  earth  noise  recorded  by  the  verti- 
cal seismographs  rises  at  a rate  of  12  to  14  db/octave  and  despite  the 
data  from  numerous  power  spectra  and  cross-spectra  and  the  information 
obtained  from  calculations  with  the  Sorrells  model,  still  only  general 
comments  can  be  made  about  the  source  of  earth  noise  recorded  by  verti- 
cal sei smographs . First  if  the  seismograph  vault  is  located  on  com- 
petent bedrock,  the  vertical  earth  noise  is  independent  of  ordinary 
atmospheric  pressure  variations.  Second  as  reported  by  Savino  e_t — a_l_. 
(1972)  the  window  in  the  earth  noise  spectra  near  40  sec  may  be  a 
"transition  between  the  swell-generated  microseisms  at  periods  less 
than  30  sec  and  non-propagating  seismic  ground  motion  of  atmospheric 
origin  at  periods  greater  than  40  sec".  A recent  paper  by  Sorrells 
and  Douze  (1974)  provides  additional  insight  into  a possible  atmos- 
pheric source  for  the  long-period  earth  noise  recorded  by  a vertical 
seismograph.  They  report  on  an  experiment  conducted  near  Dallas,  lexas, 
using  a microbarograph  array  and  s i ngl e .verti cal  seismograph  and  micro- 
barograph pair  approximately  65  km  from  the  array.  On  23  Feb.  1972 
the  microbarograph  array  detected  a very  low  level  signal  with  oscil- 
lations of  20-100  sec  period  from  the  northwest  travelling  with  an 
acoustic  velocity  of  about  330  m/sec.  The  signal  was  probably  a 
Mountain  Associated  Wave  (MAW)  train  (Larson  et_al  . , 1971)  generated 
in  the  coastal  ranges  of  British  Columbia.  For  tTTe  time  period  during 
which  the  acoustic  signal  was  being  recorded  by  the  microbarograph 
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Figure  7*  A set  of  composite  earth  noise  spectra  for  the  vertical  seismographs 


at  10  HGLP  stations,  a and  b.  Note  ZLP  - Z is  repeated  as  a reference. 
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array,  the  signal  from  the  vertical  seismograph  was  found  to  cor- 
relate with  the  signal  from  the  single  microbarograph  adjacent  to 
the  seismograph.  A very  important  factor  in  this  observation  is 
that  a condition  of  near  zero-velocity  winds  existed  in  the  vicinity 
at  the  microbarograph;  the  mi crobarometr i c noise  generated  by  the 
local  wind  would  have  masked  the  MAW-type  acoustic  waves  recorded 
by  the  single  microbarograph.  Unfortunately,  this  single  case  does 
not  provide  sufficient  data  to  establish  MAl-l  waves  as  a constant 
source  of  the  observed  long-period  earth  noise  recorded  by  vertical 
seismographs,  but  it  does  provide  additional  impetus  to  continue  a 
study  at  OGD  of  the  atmosphere  as  a possible  source  of  long-period 
earth  noise. 

Herron  et  al.  (1969)  studied  the  mi crobarographi c noise  spectrum 
using  a large  array  (250  km)  with  OGD  within  the  array.  The  pressure 
spectra  rose  at  about  12  db/octave  toward  the  long-period  end  begin- 
ning at  about  150  sec.  The  pressure  spectra  displayed  the  same  sea- 
sonal variation  that  Savino  et  al . (1972)  observed  in  the  earth  noise 
spectra  from  OGD.  These  two  observations  provide  an  indication  that 
the  two  noise  spectra,  pressure  and  earth,  are  related  but  not  neces- 
sarily as  cause  and  effect. 

The  most  helpful  data  in  determining  the  source  of  earth  noise  j 

recorded  by  horizontal  seismographs  at  the  deep  HGLP  stations  and  at 
the  shallow  and  intermediate-depth  stations  during  quiet  intervals 
are  the  shape  and  absolute  level  of  the  power  spectra  and  the  near- 
equality of  the  noise  level  observed  on  both  vertical  and  horizontal 
se i smographs . 

If  power  spectra  of  the  earth  noise  recorded  by  the  horizontal  I 

seismographs  are  compared  with  power  spectra  of  the  earth  noise  re-  | 

corded  by  the  vertical  seismograph  for  ZLP  (Figure  2),  it  is  apparent  1 

that  during  the  quietest  intervals  the  shapes  and  absolute  levels  of 
the  power  spectra  are  quite  similar.  The- same  observation  can  be  made 
for  the  other  deep  stations.  This  may  be  taken  as  an  indication  that 
earth  noise  recorded  by  the  horizontal  and  the  vertical  seismographs 
are  generated  by  the  same  source.  Cross-spectra  have  been  calculated 
between  the  earth  noise  recorded  by  the  vertical  and  that  recorded 
by  the  horizontal  seismographs  and  the  correlation  was  not  significant 
at  the  952  confidence  limits.  Thus,  if  there  is  a single  source  for 
the  two  noise  types,  it  must  be  "random"  in  effect  or  produce  only  very 
low  level  coherence  among  the  three  components. 

Finally,  the  12  to  16  db/octave  rise  in  the  earth  noise  spectra 
is  a very’  strong  indicator  that  the  principal  source  of  earth  noise 
recorded  on  both  vertical  and  horizontal  seismographs  is  not  earthquake 
coda.  Two  opposing  theories  on  earthquake  radiation  (Archanbeau , I960 
and  Brune,  1970)  state  that  earthquake  spectra  are  either  falling  off 
tov/ard  the  long-period  end  or  are  constant  at  long-periods,  but  not 
rising  at  12  to  16  db/octave.  Thus,  earthquake  coda  are  probably 
eliminated  as  a principal  source  of  earth  noise  recorded  on  vertical 
sei smographs . 

Conclusions.  The  recent  installation  of  11  high-gain  seismograpn  -* 

stations  on  a world-wide  basis  and  the  variety  of  geologic  and  tectonic  : 

settings  of  these  stations  provide  an  important  new  data  base  on  the  I 

nature  of  long-period  earth  noise.  Examination  of  several  years  of  | 


seismograms  and  analysis  of  digital  data  recorded  at  these  sites 
lead  to  the  following  conclusions  concerning  the  spectral  character- 
istics and  temporal  behavior  of  long-period  earth  noise: 

1.  In  general,  earth  noise  spectra  calculated  for  all  the  sites 
exhibit  a minimum  level  in  the  period  range  25  to  45  seconds. 

2.  The  level  of  earth  noise  recorded  by  vertical  seismographs 
is  apparently  independent  of  local  meteorological  conditions.  The 

low  noise  levels  achieved  at  the  shallow  sites  are  a result  of  the  sites 
being  located  on  competent  rock  formations  and  having  rigid  environ- 
mental control. 


3.  The  level  of  earth  noise  recorded  by  horizontal  seismographs 
is  dependent  upon  the  depth  of  overburden.  At  the  shallow  sites 
(depth  of  overburden  < 1 00  in)  the  earth  noise  levels  exhibit  a general 
diurnal  cycle.  This  cycle  appears  to  be  caused  by  local  tilting  of 
the  ground  associated  with  daily  variations  of  atmospheric  turbulence. 


4.  The  special  UGLP-type  horizontal  seismographs  on  the  103  and 
153  rn  level  of  the  Ogdensburg  Deep  Mine  Observatory  and  the  other  HGLP 
horizontal  seismographs  show  that  the  earth  noise  (tilts)  associated 
with  atmospheric  turbulence  are  attenuated,  in  general,  to  about  10/6 
of  the  surface  level  by  about  150  m in  competent  rock  formations, 
i.e.,  noise  level  at  150  m depth  is  10%  of  the  noise  level  on  the  surface 

5.  The  source  of  the  long-period  earth  noise  recorded  by  the  verti- 
cal s e i smog r a phs  and  by  the  horizontal  seismographs  at  the  deep  station 
is  probably  of  atmospheric  origin.  This  aspect  of  the  problem  of  earth 
noise  source  will  receive  further  investigation. 


2 . Attenuation  of  Long-Period  Earth  Noise  with  Depth 

A vertical  array  of  long-period  seismographs  was  constructed  in  the 
Ogdensburg  Deep  Mine  Observatory  to  measure  the  attenuation  of  earth 
noise  with  periods  near  40  sec  over  the  depth  range  0 to  543  m.  Using 
a model  derived  by  Sorrells  (1971)  to  predict  ground  motions  produced 
by  a slow -moving,  plane-wave  pressure  disturbance,  we  calculated  the 
expected  attenuation  with  depth.  According  to  the  model  this  attenua- 
tion is  dependent  upon  the  velocity  at  which  the  pressure  disturbance 
is  propagating.  The  observed  atteivuations  correspond  to  velocities 
that  are  consistent  with  the  velocities  of  atmospheric  phenomena  tnat 
have  previously  been  reported  to  produce  ground  motions  recordec  by 
long-period  seismographs. 

The  tilt  motions  recorded  by  horizontal  seismographs  appear  to  be 
generated  by  low  velocity  (<5  to  10  rn/s)  turbulence  in  the  atmospheric 
boundary  layer  immediately  above  the  surface  of  the  earth.  These 
motions  are  significantly  attenuated  by  about  100  to  150  m of  compe- 
tent overburden  (Figure  0).  The  vertical  and  horizontal  displacements 
and  tilt  motions  produced  by  atmospheric  disturbances  travelling  at 
acoustic  velocities  (about  300  m/s)  are  not  significantly  attenuated 
by  500  m of  overburden.  Thus,  the  signal -to-noise  ratio  on  horizontal 
seismographs  can  be  significantly  improved  by  burial  at  depth  of  several 


hundred  meters  which  removes  the  major  portion  of  the  tilt  noise 
associated  with  low  velocity  winds,  but  the  signal-to-noise  on  the 
vertical  and  horizontal  seismographs  cannot  be  further  improved 
without  going  to  impractical  depths  of  burial  (more  than  two  or 
three  km). 


3 . Detection  of  Air  Waves 

Acoustic-gravity  waves  from  an  event  on  14  October  1970,  pre- 
sumed to  be  a Chinese  atmospheric  explosion,  were  observed  on  high- 
gain  vertical  and  horizontal  component  seismograms  written  at  three 
different  locations  in  the  world  (Figure  9).  Reliable  group  velocity 
data  for  the  atmosphere  over  propagation  paths  as  long  as  97°  (10,800 
km)  were  obtained  for  the  period  range  of  30  seconds  to  375  seconds. 
These  data  are  in  agreement  with  theoretical  dispersion  curves  that 
are  based  on  the  COSPAR-mode 1 atmosphere  with  the  effects  of  winds 
included.  Seismic  body  and  surface  waves  from  this  atmospheric  event 
were  also  recorded  at  the  three  stations  and  were  used  to  determine 
an  epicenter. 

There  are  two  important  points  about  the  recordings  of  the  acoustic 
gravity  waves  on  the  seismograms  in  Figure  9.  First,  because  of  the 
very  long-period  character  of  these  waves,  it  should  be  quite  simple 
to  dptpct  thpir  presence  In  the  midst  of  shorter-period  surface  waves 
from  a seismic  event  by  applying  low-pass  filtering.  Second,  these 
waves  are  recorded  with  a very  high  signal-to-noise  ratio  in  the  peri.od 
range  60-100  sec.  A part  of  the  seismic  background  is  included  on 
each  seismogram  for  comparison.  These  acoustic  gravity  waves  were 
also  recorded  on  mi crobarograms  at  stations  near  Palisades,  New  York, 
approximately  50  km  SSE  of  Ogdensburg  and  at  College,  Alaska.  Com- 
parison of  the  mi crobarograms  with  the  high-gain  seismograms  from  0GD 
and  FBK  shows  that  the  infrasonic  waves  with  periods  between  60  and 
100  sec  were  recorded  w'th  higher  signal-to-noise  ratios  on  the  seis- 
mograms than  on  the  mi crobarograms . At  periods  of  >100  sec  the  response 
of  the  high-gain  seismographs  falls  off  rapidly  with  period  (^18  d b / 
octave),  and  thus  the  recorded  amplitude  of  the  infrasonic  waves  is 
reduced  at  these  longer  periods. 

On  14  April  1971,  slow  speed  gravity  waves  from  another  pressure 
disturbance  of  meteorological  origin  were  observed  with  periods  between 
240  seconds  and  360  seconds  on  seismograms  at  Ogdensburg,  Mew  Jersey 
(Figure  10).  The  seismic  recordings  of  the  acoustic-gravity  and  the 
slow  speed  gravity  waves  are  attributed  to  ground  motion  produced  by 
surface  loading  and  not  to  direct  pressure  effects  on  the  instruments. 
The  observed  displacements  and  tilts  from  both  events  are  in  agree- 
ment with  those  predicted  by  static  loading  theory.  At  the  Ogdensburg 
station,  the  rock  layer  above  the  instruments  (543  m)  acts  as  a wave- 
length filter  to  suppress  wind  noist  of  short  wavelength  and  to  ^.barica 
signals  from  long  wave-length  (coherent)  disturbances  in  the  atmos- 
phere . 
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Fig.  Seismograms  showing  long-period  body  phases,  surface  waves,  and  acoustic  gravity 
waves  for  the  presumed  Chinese  atmospheric  explosion  of  October  14,  1970,  from  the  vertical 
(ZHI)  and  horizontal  ( N III,  EHI)  component  high-gain  seismographs  at  Ogdcnsburg  (A  = 
10,798  km),  Charters  Towers  (A  = S939  km),  and  Fairbanks  (A  = 7276  kin)  and  from  the 
long-period  vertical  WW’SSN  component  seismometer  at  Ogdcnsburg.  The  time  base  on  all 
seismograms  is  1 hr/line. 
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Fig.  10.  Mirrobarogrnm  and  three-component  high-gain  and  WYVSSN  seismograms  for  a local 
atmospheric  event  on  April  11,  1971,  recorded  in  the  mine  observatory  at  OGD.  The  scale' 
given  for  the  microbarogram  refers  to  the  sensitivity  of  the  microbarograph  at  periods  be- 
tween about  GO  and  2S0  see.  On  either  side  of  this  period  range  the  response,  of  the  miero- 
barograidi  falls  olT  at  about  6 db/octave. 


4 . Very-Lonn -Period  Mi crosei sms 

The  double  frequency  (OF \n\c\ose] \sml 0nPi  i near°i nteracti on°of 
sec)  are  surface  waves  generated  by  the  non-  *|]et"J®ggi ns , 1950; 

trains  of  oppositely  travelling  o (pri  mi  crosei  sms  (period 

Hasselmann,  1963).  Jhe  primary  ■ fr^  ^ by  llaubrich  and  IlcCamy 
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for5 transient0 background  noise  P'-nonena  are^n^as^eU^nder^ 

^ic.,°^rSca^n.u  ddp  srace -cs 4--ated  tb0dwaken0f 

an  unusually  swift  tropical  storm  far  at  seaman  s,ow.5peed 

(1972)  studied  the  srou"d  ^placen^  p J metc ecological  origin, 
gravity  waves  from  a pressure  “^U-La-  - 

Recently,  Murphy  and 

of  35  to  40  sec  Rayleigh  wave  mi croseism s Jhe  Kon”s5erg  station 

high-gain,  long-period  seismograph  * 4 Jrd  DF  and9 PF  microseisms 
is  well  situated  in  s°uthern  Norw  y ^£0  the  International 

generated  by  gales  in  the  gojth  A 2,  tiie  p0V/er  spectra  was  Sal- 

Seismic  Month  (20  Feb.  - 19  March,  1972),  the  P seVRre 

culated  for  nearly  every  day.  For  the  days^u  9pa  disp1ay  a 
DF  and  PF  microseisms  were  absent,  * P™  sec"  but  for  four  days 
pronounced  minimum  between  croseism  storms  were  present  the 

E^is2-“ii£hi::r”o55-ssJ-p«3'ssr?i::i,5..  raised  about » 

The  coherence  between  the  v^5^“\p"5n3°i;flYnic?oseismSstorms 

nnsim-t  "tie1  id  ? .in:.thc,.5Ps.rtr-;ne2’n.t 

although  the  phase i spectrum  between  the  same^  expect;,d  for  Rayleigh 
show  the  90°  or  270  relation  bet<  ™rd  9„.  „!!  the  period  decreases 
waves,  there  is  a trend  in  the  a does  not  ruie 

The  departure  of  the  i°ast  two  reasons.  The  phase 

out  Rayleigh  wave  type  grou 01  possible  changing  phase  shift 

ini  mott on°i s PrSlf U^rTo-bln./jIth  significant  Love  wave 
motion,  uni ty  coherence  is  not  expected.) 

' The  signal  from  the  three filter™ 
sample  was  passed  through  1 t 0f  these  traces  shows  three 

with  a corner  period  of  3b  sec.  ...  P j,Ptween  the  vertical  and  M/S 
things:  rc a s °n a t> 1 y g °od  c ° r r e 1 a tl  vertical  trace  lead 

^^Hh^ertical  - */«  traces  correlate 

far  less  well. 

There  was  no  significant  coherence  be t«eenT Jjeg/three  observa- 
tions^?0"^ gnifica^  component 

wi th^verti calCl eadi ng^  and  no  significant  coherence  between  vertica 
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and  E/W  components,  are  good  indicators  that  what  we  have  observed 
are  35  to  40  sec  Rayleigh  wave  mi crosei sms  with  retrograde  motion 
coining  from  a source  north  of  KON.  A check  of  the  surface  weather 
charts  for  the  northern  hemisphere  shows  that  there  is  a large 
low  pressure  system  north  of  KGU  for  the  proper  time  interval,  t It  a t 
is,  allowing  for  ocean  wave  propagation  from  the  storm  to  the  coast 
line.  It  is  likely  that  this  storm  is  the  source  of  the  ocean  wave 
energy  responsible  for  the  t)F  and  PF  microseism  storm.  The  pos- 
sibility that  the  observed  signal  was  the  coda  of  a large  earthquake 
was  eliminated  by  checking  the  list  of  earthquakes  for  the  Inter- 
national Seismic  Month  and  by  searching  for  similar  noise  at  the 
same  time  at  the  Toledo  and  Ogdensburg  HGLP  stations.  No  large 
earthquake  is  listed  near  the  time  of  the  mi  crosei  sin  storms,  and 
the  signals  were  not  observed  as  far  away  as  TuL  or  *'60. 


MAGNITUDE  THRESHOLDS  OF  THE  IIGLP  STATIONS 

1 . Detection  Capabilities  of  8 High-Gain  Stations  versus  the  WWSSN 


A detailed  investigation  of  seismograms  from  eight  high-gain 
stations  for  the  time  period  8 September  1971  to  15  September  1971 
(inclusive)  was  undertaken  to  compare  the  number  of  earthquakes 
recorded  at  these  eight  stations  with  the  number  reported  by  NOAA 
in  the  monthly  listings  of  the  PDE.  While  long-period  surface  waves,' 
of  both  the  Rayleigh  and  Love  type,  are  the  primary  type  of  seismic 
wave  that  the  high-gain  detection  results  are  based  on,  events  are 
reported  in  the  PDE  listings  on  the  basis  of  the  recording  of  short- 
period  body  phases  at  at  least  a few  of  the  stations  of  the  WWSSN. 

The  results  are  shown  in  Figure  11  in  the  form  of  a histogram. 

The  numbers  inside  the  hatched  (lower)  sections  in  Figure  11 
refer  to  the  number  of  earthquakes  reported  in  the  PDE  and  observed 
at  the  station  indicated;  the  solid  sections  refer  to  the  number 
masked  at  a particular  station;  and  the  stipled  sections  to  the 
number  of  unreported  events  observed  at  a particular  high-gain 
station.  Of  the  103  reported  events  (the  hatched  section  of  the 
column  designated  WWSSN),  79  were  recorded  on  one  or  more  of  the 
seismograms  from  the  eight  stations  as  indicated  in  the  column  adja- 
cent to  that  for  the  WWSSN.  Of  the  remaining  24  reported  earthquakes, 
15  were  definitely  not  observed  at  any  of  the  high-gain  stations. 

Six  of  these  15  events  were  assigned  focal  depths  greater  than  100  km, 
6 others  were  not  assigned  any  body-wave  magnitudes,  and  the  remaining 
3 were  assigned  magnitudes  less  than  4.0  and  occurred  at  distances 
greater  than  25  to  30°  from  the  closest  high-gain  station. 

The  number  of  events  observed  on  the  combined  high-gain  seismo- 
grams but  unreported  by  PDE  is  73.  This  gives  a total  of  152  events 
observed  during  this  time  period  or  approximately  20  per  day.  During 
this  same  time  period  short-period  body  phases  from  197  events  were 
reported  in  the  bulletin  from  the  Large  Aperture  Seismic  Array  (LASA) 
in  Montana,  ilote  that  we  cannot  compare  these  numbers  with  one  from 
the  PDE  program  since  many  of  the  earthquakes  that  are  detected  at 


a few  isolated  stations  of  the  WWSSN  are  not  located  and  thus  go 
unreported.  For  this  reason  the  number  of  unreported,  but  de- 
tected, events  for  the  WWSSN  is  left  open-ended.  To  be  sure  it 
is  in  excess  of  the  total  number  observed  on  the  high-gain  seismo- 
grams and  probably  also  the  number  reported  in  the  LASA  bulletin. 

The  number  of  NOAA-repor ted  earthquakes  with  surface  wave 
trains  masked  by  the  coda  of  larger  events  is  9 for  the  combined 
seismograms  from  all  eight  high-gain  stations.  This  number  is 
not  greatly  different  from  the  numbers  of  masked  events  at  any 
one  station  and  points  up  one  of  the  most  important  problems  with 
the  detection  of  long-period  surface  waves  from  shallow  events. 

While  more  sophisticated  analysis  techniques  than  the  visual  one 
employed  in  this  study  would  probably  identify  some  of  these 
masked  events,  it  remains  to  be  seen  how  successful  these  techniques 
will  be  in  routine  practice.  Notice,  however,  that  the  number  of 
surface  w aves  (79  + 73  = 153)  detected  on  the  combined  high-gain 
seismograms  is  about  double  the  number  detected  by  a single  high- 
gain  station. 

2 . Regional  Magnitude  Threshold 

Seismograms  from  the  stations  at  OGD,  CTA,  FBK,  EIL  and  CHG 
were  analysed  for  more  than  1000  shallow  focus  (h^CO  km)  earth- 
quakes reported  in  the  PDE  to  assess  the  dependence  of  the  detection  , 
capability  for  long-period  surface  waves  of  a particular  station 
on  the  geographical  locations  and  epicentral  distances  of  earth- 
quakes. Figure  12  shows  the  surface  wave  detection  thresholds  with 
respect  to  a particular  station  for  59  seismically  active  areas. 

The  thresholds  are  expressed  in  terms  of  body-wave  magnitude, 
m,  , as  determined  by  NOAA  (PDE)  with  the  exception  of  western  North 
America  (exclusive  of  Alaska)  where  Evernden's  (1967)  formulas  for 
mh  were  used.  Magnitudes  of  earthquakes  outside  western  North 
America  are  not  likely  to  change  significantly  if  Evernden's  formu- 
las are  used  rather  than  those  of  NOAA  since  the  percentage  of  sta- 
tions at  distances  less  than  18°  is  nearly  always  small.  The  values 
of  nil  indicated  on  the  map  represent  the  smallest  magnitude  earthquakes 
frorn  which  surface  waves  were  consistently  observed  at  one  of  the 
high-gain  stations  during  the  particular  time-period  of  this  study. 

When  a value  of  the  magnitude  threshold  is  given  in  Figure  12  for  a 
particular  region,  it  means  that  all  earthquakes  reported  in  the  PDE 
(the  data  base)  for  that  region  with  magnitudes  equal  to  or  greater 
than  this  threshold  value  generated  observable  surface  waves. 

The  less  than  symbol  (<)  indicates  that  we  were  only  able  to 
determine  an  upper  limit  for  the  value  of  the  magnitude  threshold 
in  a particular  region  because  of  the  inadequate  number  of  small 
magnitude  earthquakes  reported  in  the  PDE  for  that  region.  Improved 
detection  and  reporting  of  short-period  body  phases  would  allow  the 
determination  of  more  accurate  surface-wave  thresholds  for  events  in 
these  regions. 
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The  regions  determined  using  OGD  data  include  the  North  and 
South  Atlantic  Ocean,  the  Caribbean,  South  America,  western  United 
States,  and  Vancouver  Island.  Six  months  of  data  from  F B K were 
used  to  determine  thresholds  in  Alaska,  the  Aleutians,  and  the 
Kuri  1 es-Kamchatka  regions;  four  months  of  C I i G data  for  Hokkaido, 
Honshu,  Ryukyu  and  the  mid-Indian  Rise;  six  months  of  LIL  data 
for  Tibet,  Iran,  and  Turkey;  and  six  months  of  CTA  data  for  the 
remaining  regions  of  the  western  Pacific  and  Easter  Island.  It 
is  important  to  note  that  in  many  of  these  regions,  surface  waves 
from  earth quakes  with  mb  values  as  much  as  0.75  units  less  than  the 
threshold  values  indicated  in  Figure  12  were  observed. 

The  better  determined  threshold  values  given  in  Figure  12  are 
plotted  in  Figure  13  as  a function  of  epicentral  distance.  In  the 
distance  range  20°  to  25°  , the  detection  threshold  for  long- 
period  surface  waves  is  approximately  nv  = 4.1  (PDF.).  The  dashed 
curves  in  Figure  13  represent  the  expected  thresholds  under  relatively 
high  and  low' noise  conditions.  These  curves  were  computed  using  the 
magnitude  relationships  described  by  Gutenberg  (1945).  The  noise 
level  of  15  nig  represents  a record  amplitude  of  1 mm  at  20  s which  is 
about  the  lowest  limit  for  visual  analysis  under  ideal  conditions. 
While  these  thresholds  for  surface  waves  with  periods  of  20-40  s 
represent  an  improvement  over  existing  seismograph  systems,  the 
various  digital  enhancement  techniques  that  can  be  applied  to  the 
data  from  the  high-gain  systems  should  result  in  a further  lowering 
of  the  threshold  values. 


3 . Detection  Capabilities  of  the  HGL.P  Stati  ons  During  International 
Seismic  Month 


Recently  the  Seismic  Discrimination  Group  at  Lincoln  Laboratory 
compiled  a special  listing  of  the  occurrence  of  as  many  events  as 
possible  for  the  time  period  20  February  1972  through  19  March  19/2. 
Hereafter  this  special  listing  and  the  above  1 month  time  period 
will  be  referred  to  as  the  International  Seismic  Month. (ISM).  The 
final  version  of  the  ISM  listing  contains  996  entries  identified  on 
the  basis  of  short-period  body  waves  recorded  on  the  various  arrays 
(LASA,  N0RSAR,  ALPA);  the  Canadian,  Swedish,  and  United  Kingdom 
seismograph  networks;  and  the  WWSSN  stations.  Tl  i;  ^ ISM  1 i sting  serves 
as  an  excellent  data  base  to  be  used  in  the  determination  of  approxi- 
mate detection  thresholds  of  the  high-gain  stations.  The  results, 
discussed  in  this  study  were  based  on  a preliminary  listing  conta-in- 
ing  840  events. 


Figure  14  shows  the  percentage  of  reported  events  above  the 
plotted  magnitude  for  which  surface  waves  were  observed  at  one  or 
more  of  the  high-gain  stations.  Only  shallow  (h<100  km)  events, 
which  had  two  or  more  stations  reporting  a magnitude  were  used  in 
this  comparison.  Note  that  surface  waves  were  observed  for  90% 
of  the  reported  events  of  m.  4.6  or  greater.  Half  of  the  1 0 a of 
these  events  not  observed  were  the  result  of  masking  by. the  coda 
of  earlier  larger  events.  Significantly,  the  largest  single  event 
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not  in  coda  for  which  surface  waves  were  not  observed  was  of 
4.8,  and  one  additional  event  of  m.  4.7  was  also  missed.  All 
the  remaining  events  not  in  coda  for  which  surface  waves  were  not 
observed  were  of  m^  4.6  or  less. 

Two  other  important  observations  that  are  not  obvious  from 
Figure  14  come  from  the  analysis  involved  in  the  preparation  of 
this  figure.  They  are  the  following: 


1.  Approximately  10%  of  the  events  observed  at  the  deep 
stations  E1L , OGD,  and  KON  were  observed  on  the  basis  of  long- 
period  Love  waves  only.  The  reason  for  restricting  this  state- 
merit  to  the  deep  stations  is  that  the  horizontal  component  seis- 
mographs  at  these  sites  operate  with  peak  magnifications  equi- 
valent to  those  of  the  vertical  component  systems.  In  addition 
the  low  noise  levels  on  the  horizontal  component  seismograms 
from  these  three  stations  do  not  undergo  any  obvious  diurnal 
variations  as  do  the  horizontal  noise  levels  at  the  other 
shallower  sites  (Murphy  et  al . , 1 972  ; Savino  et  a 1 . , 1 972a  and  b). 


2.  Surface  waves  with  periods  between  30  and  50  seconds  were 
recorded  from  more  than  90%  of  the  465  observed  events.  This 
percentage  is  greater  than  that  for  the  number  of  events  for  which 
surface  waves  with  periods  near  20  seconds  we re  observe  . 


The  histogram  in  Figure  15a  gives  an  idea  of  how  the  VLP 
stations  as  single  stations  (e.g.,  CTA  column)  fared  in  t^nis  of 
the  ISM  listing  (column  headed  ISM).  The  numbers  inside  the  lightly 
stippled  sections  of  each  of  the  VLP  columns  refer  to  the  number  of 
world-wide  earthquakes  reported  in  the  ISM  that  excited  surface 
waves  of  either  the  Rayleigh  or  Love  type  with  periods  between  about 
15  to  60  seconds  and  which  were  observed  at  that^p a rticular  (see 
column  heading)  VLP  station.  Mote  that  about  jo*  (J6j)  of  the  840 
reported  events  were  detected  at  at  least  one  of  the  VLP  stations, 
more  than  twice  as  many  than  the  number  detected  at  any  one  site. 
While  it  is  not  entirely  evident  from  this  figure,  the  combination 
of  the  five  stations  EIL,  FBK , CTA,  ALQ , and  KON  are  responsible  for 
more  than  90%  of  the  465  observed  events. 

The  varying  heights  of  the  individual  VLP  columns  in  ^ipur<^ . 
reflect  the  amount  of  time  no  records  were  available  from  the  differ' 
ent  stations.  For  all  the  stations  exceptKIP  and  CHG , these  times 
were  mainly  due  to  analog  record  changing  intervals;  digital  data 
are  available  for  these  times.  The  KIP  station  became  operational 
on  7 March  1972,  more  than  half  way  through  the  ISM  time  period. 

This  accounts  for  the  relatively  small  number  (*: 01  ) of  events 
considered.  Approximately  7 days  of  recordings  from  CHG  were  no 

available. 

The  numbers  in  the  open  sections  of  the  VIP  columns  and  the 
additional  numbers  in  the  hatched  sections  of  the  JLO,  OGD,  and 
KON  columns  designate  the  numbers  of  reported  events  that  we  - 
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TABLE  2. 


90%  surface  wave  detection  thresholds  (rnu) 

b" 


20°  to  30° 

30°  to  40° 

80°  to  100( 

CTA 

4.0 

4.9 

EIL 

4.2 

4.9 

ALQ 

' 4.2 

f — 

4.4 

KON 

4.0 

• 4.5 

TLO 

^ 4.5 
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masked  by  either  body  waves  and/or  surface  waves  (coda)  from  a 
larger  event  (open  sections)  or  masked  by  large  amplitude  micro- 
seism  storms  comprised  of  surface  waves  with  periods  near  H 
20  seconds.  Note  that  the  severity  of  the  masking  problem  ranges 
from  46%  of  the  events  considered  at  KON  to  16.4%  at  FBK.  The 
effect  of  combining  all  nine  VLP  stations  is  to  reduce  the  number 
of  masked  events  to  100  out  of  340  or  12%.  Actually  the  number 
or  percentage  of  masked  events  based  on  all  the  VLP  stations 
(12%  here)  for  the  most  part  will  not  be  significantly  smaller 
than  the  smallest  number  or  percentage  obtained  at  some  one  sta- 
tion (e.g.,  FBK  137).  This  is  a result  of  the  well  developed, 
and  widely  recorded,  coda  of  surface  waves  associated  with  shallow 
events  of  magnitudes  greater  than  5.0  to  5.5.  During  the  ISIT 
month  two  events  with  magnitudes  greater  than  6.0  occurred  and 
were  responsible  for  50  of  the  100  events  masked. 


During  the  winter  months  in  the  northern  hemisphere,  and  es- 
pecially in  the  north  Atlantic,  microseism  storms  of  both  the  double 
frequency  (7  to  9 sec)  and  primary  frequency  type  (15  to  18  sec)  are 
particularly  common.  For  instance  microseism  storms  that  occurred 
during  the  ISM  month  accounted  for  as  much  as  33%  of  the  visual  mask- 
ing problem  at  KON.  In  this  regard,  however,  it  is  interesting  to 
note,  although  well  understood  (Haubrich  and  McCamy,  1969  ; liassel- 
m ann,  1963),  that  while  several  double  frequency  storms  were  observed 
on  the  low-gain  or  standard  seismograms  from  KIP,  because  of  the 
structure  of  the  Hawaiian  Islands  (the  absence  of  a sloping  shelf) 
there  were  no  observable  primary  frequency  storms.  Thus,  in  terms 
of  masking,  the  long-period  (T>15  sec)  detection  thresholds  deter- 
mined during  the  ISM  at  this  mid-ocean  station  depended  only  on  the 
coda  problem  discussed  above. 


As  both  nicroseism  storms  and  the  occurrence  of  large  magnitude 
earthquakes  are  transitory  phenomena  it  is  of  interest  to  compare 
the  performance  of  the  nine  VLP  stations  during  a time  interval 
when  both  of  these  masking  sources  are  active  with  a time  interval 
when  both  sources  are  absent.  In  Figure  15b,  the  performance  of  the 
VLP  stations  during  two  48  hour  time  periods  within  the  ISM  month 
are  compared.  The  time  interval  on  the  left-hand  side  corresponds 
to  the  absence  of  large  earthquakes  and  microseisms  while  that  on 
the  right-hand  side  corresponds  to  the  dominance  of  one  or  both  sources. 
Note  the  drastic  change  in  the  detection  capabilities  from  72  outof  J 
(80%)  during  quiet  times  to  19  out  of  90  (21%)  reported  events  being  ob- 
served when  masking  is  severe.  The  corresponding  percentage  change  in  tne 
number  of  masked  events  is  tenfold.  While  masking  is  obviously  a 
problem  with  the  VLP  systems,  it  must  again  be  pointeooutt natal 
the  results  presented  so  far  are  based  only  on  visual  analysis  ot 
nnalnu  recordings.  It  is  to  be  expected  that  some  rather  simple, 
filtering  techniques  (low-pass,  band-pass)  and  the  more  sophisticated 
tpchn^oues^ofmatchedfilterinq  and  polarization  filtering  (Choy  and 
1973  that  can  be  applied  to  the  digitally  recorded  data  at 
HcCany,  J improve  (lower)  the  detection  thresholds,  especially 

S^g'ecord^  1 n^a,  s ?ha  t are  visually  dominated  by  either  micro- 
seisms  or  interfering  events. 


A final  point  to  be  considered  concerning  the  results  in  Figure 
15  the  question  as  to  how  do  the  VLP  stations  fare  detection- 
wise  when  only  shallow  (here  taken  as  h<50  km)  events  are  considered. 
. 3 ,r?  iCJ31Vfs>t  ,e  answer  to  this  question.  Comparing  the  right- 
and  (all  depths)  and  left-hand  side  of  this  figure,  we  find  that 
the  percentages  of  events  observed  (SIG),  not  observed  (MO  SIG), 
and  masked  do  not  change  when  the  ISM  listing  is  screened  for  deeper 
events.  To  some  extent  this  result  is  surprising.  One  possible 
explanation,  however,  is  that  depths  greater  than  normal  ( h >33  km) 
were  in  many  instances  not  assigned  to  the  very  small  magnitude 
events.  This  is  not  unexpected  since  the  very  small  events  were 

probably  poorly  recorded  making  proper  depth  determination  extreme^ 
difficult. 


To  study  detection  thresholds  as  a function  of  distance,  all 
those  events  listed  in  the  ISM  as  having  a focal  depth  of  50  km  or 
ess,  and  assigned  a body-wave  magnitude  (nik) , were  considered.  In 
Table  2 values  of  the  90S  detection  thresholds  in  two  different  10° 
distance  ranges  are  given  in  terms  of  ni.  for  5 of  the  VLP  stations 
because  of  insufficient  data  at  CTA  and  EIL  in  the  20°  to  30°  distance 
range  threshold  values  were  instead  determined  in  the  30°  to  40° 
range  at  these  two  staions.  From  Table  2 we  find  that  mean  values 
°I  ^etection  thresholds  for  surface  waves  range  from  m,  =42 

at  30  to  about  m5  = 4.7  at  90°.  It  should  be  noted  again  that  these 
threshold  values  are  based  on  visual  analysis  of  seismograms  and 
that  the  filtering  techniques  described  by  Choy  and  McCamy  (1973) 
should  lower  these  values  by  at  least  0.2  mb  units. 


T N S r ? N T R A f T A ? T BETWEEN  EARTH()UAKES  AMD  PRESUMED  UNDERGROUND  EXPLOSIONS 


Previous  results  from  a high-gain  station  on  the  Mc-mu  discrim- 
inant between  earthquakes  and  underground  explosions  in  different 
regions  of  the  world  were  reported  by  Molnar  e t a 1 . ( 1 969  ),  Evern- 
den  e t a 1 . (1971)  and  Savino  e t a 1 . (1971).  These  earlier  studies 
were  based  on  data  recorded  at  OGD.  In  this  section  we  will  concen- 
trate on  results  from  one  of  the  more  recent  stations,  that  at 
Eilat,  Israel.  This  station  is  approximately  30®  from  presumed  Rus- 
sian underground  explosions  and  exhibits  a lower  detection  threshold 
for  surface  waves  from  these  presumed  explosions  and  earthquakes  in 
the  same  general  region  than  the  stations  at  Chiang  Mai,  Thailand 
(CHG);  Kongsberg,  Norway  (KON);  and  Toledo,  Spain  (TLO). 

In  Figure  16,  peak- to-peak.  amplitudes  of  Rayleigh  waves  recorded 
at  Eilat  with  periods  near  20  sec  (left-hand  plot)  and  40  sec  (right- 
hand  plot)  from  earthquakes  in  the  Ki rzi g-S i nki ang-Tadzh i k-Ti bet 
region  (closed  circles)  and  presumed  Russian  explosions  (open  circles) 
aie  plotted  as  a function  of  mb  (NOAA).  The  ordinate  on  the  20  sec 
plot  describes  the  Ms  scale  for  these  events  as  indicated  by  the 
numbers  and  arrows.  Although  the  data  are  limited  there  are  two 
rather  important  points  about  Figure  16.  Firstly,  the  earthquake 
and  explosions  populations  are  completely  separated  even  though  the 


T = 20  SEC 


T « 40  SEC 


Fig.  If.  M,~mb  at  20  and  40  s (Rayleigh  waves)  for  earthquakes  in  the  Kirgiz- 
Sinkiang-Tadzhik-Tibet  region  (closed  cireles)  and  presumed  Russian  explosions 
(open  circles).  All  of  these  measurements  were  taken  f-om  seismograms  from  the 
station  at  Eilat,  Israel.  Note  the  pronounced  separation  of  the  two  populations. 
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Fig.  17.  The  rejection  exponent  N of  the  filter 
transfer  function  determines  I he  amount  of  re- 
jection of  signals  whose  angle  0 arrives  incorrectly 
polarized  or  off  azimuth. 
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in e a s u re ino n ts  arc  restricted  to  a single  station.  Secondly,  the  two 
populations  do  not  show  any  sign  of  convergence  at  the  low  magnitudes, 
although  it  must  be  remembered  that  tl  NOAA  (PDE)  reports  do  not  ^ 
include  most  events  from  this  area  smaller  than  about  = 4.5  to  4.u. 
This  latter  result  is  similar  to  that  found  at  OGD  for  earthquakes 
and  explosions  in  the  western  United  States,  the  Aleutians,  and  the 
Novaya  Zemlya  region  (Holnar  et  a 1 . , 1 969;  Savino  e t a 1 ._ , 1971). 


D.  ENHANCEMENT  OF  LONG-PERIOD  SIGNALS  BY  TIMF-VARYING  ADAPTIVE  FILTERS 


As  detection  thresholds  are  lowered,  the  performance  of  any  net- 
work becomes  limited  not  only  by  ambient  noise  but  also  by  the  mask- 
ing of  small  events  by  larger  ones.  That  is,  if  the  effect  of  back- 
ground noise  is  reduced,  the  only  competition  against  detecting  a 
desired  signal  is  the  presence  of  other  signals.  To  help  overcome 
this  limitation  on  the  performance  of  the  high-gain  stations,  we  ap- 
plied  time-varying  adaptive  filters  to  the  seismograms  from  individual 
stations.  A s i gnal -to-noi se  enhancement  of  more  than  6 db  (0.3  mag- 
nitude unit)  was  attained  for  Rayleigh  waves  from  small  events.  Love 
ivrvtf*  may  olay  an  increasingly  important  role  in  discriminating 
between  seismic  sources.  Hence  we  also  investigated  the  effect.veness 
of  using  time-varying  adaptive  filters  for  extracting  Love  waves 
masked  by  the  coda  of  a larger  events. 


1.  Filters 


The  filter  transfer  functions  are  shaped  to  pass  frequency  bands 
where  the  s i gnal -to-noi  se  ratio  is  momentarily  high,  lie  define  a 
high  s i gnal -to-noi se  ratio  in  two  different  ways:  by  the  apparent 

polarization  of  particle  motion  and  by  the  apparent  azimuth  n arriva 
The  first  criterion  is  used  in  the  'polarization  filter.  The  re  - 
tive  phase  angle  between  the  vertical  and  horizontal  motion  o(f)  is 
lispoSted  .s  1 function  of  frequency,  and  then  the  amp  itude  response 
of  the  filter  is  shaped  to  pass  either  body  waves  (o(f)  0 or  tt)  or 

Rayleigh  waves  (e(f) • +n/2).  The  second  criterion  is  used  in  the 
' 37iniuthal  ' filter.  Here  tne  apparent  direction  ul  arrival  o.  L - 
waves  on  the  two  horizontal  instruments  is  inspected  as  a function  of 
frequency.  The  amplitude  response  of  the  filter  is  shaped  to  pass 
energy  in  the  frequency  bands  where  the  apparent  signal  is  near  the 
azimuth  of  the  suspected  event. 


He  briefly  outline  the  strategy  of  our  computations.  We  are  given 
two  digitized  seismograms  x.j  and  y^  (i  =0,  1,  2,  ....  in  secon(s). 

For  the  polarization  filter,  x and]y  are  the  vertical  and  horizontal 
seismograms.  For  the  azimuthal  filter  they  are  the  two  horizontal 
seismograms.  IJe  take  segments  of  length  m and  operate  on  them  in 
turn.  In  the  kth  segment  we  take 


*ik  = x(i+km/2) 


i = 0,  1 m - 1 


yik  = y (i+km/2) 


1 


We  compute  their  Fourier  transforms  $j  and  Yj . For  the  polarization 
filter  the  phase  angle  on-  is  the  argument  of  the  cross  spectrum 

6j  = arg  K-jV-j* 

For  the  azimuthal  filter  the  angle 

Oj  = a - tan”1 (Yj/Xj ) 

is  the  difference  between  the  expected  azimuth  of  arrival  of  the 
suspected  event*  ct,  and  the  apparent  azimuth  of  arrival.  The  filter 
transfer.,  f uncti  on  for  the  kth  segment  is  then  generated: 

Fj  = sin'  o-j  for  the  Rayleigh  polarization  filter  and  the  azimuth 
filter  or  cos  o -j  for  the  body  wave  filter.  The  exponent  II  determin< 
the  amount  of  rejection  of  the  filter  for  frequencies  whose  phase 
angles  are  either  incorrectly  polarized  or  that  arrive  off  azimuth. 
Figure  17  shows  how  power  is  rejected  as  a function  of  0 for  the 
different  values  of  II.  The  Fourier  transforms-  are  then  multiplied 
by  this  filter  to  give  the  transforms  of  the  filtered  segments: 


V - XiFj 


Yj'  = YjFj 


The  first  half  of  their  inverse  Fourier  transforms,  x j'00  and  y -j  1 ( ^ , 
is  combined  with  the  last  half  of  the  preceding  (k  - 1)  segment  to 
form  the  tine-varying  filtered  output: 


*i“  - vr(k)  + <’  - wi)x(1+w/2)'0-.-i) 
y."  = wiyp<k>  + (l  - w1  )yt1+m/2)' (k-1) 


i = 0 , 1 , . . . , m/2 
U.j  = 2i/m 


Our  experience  lias  been  that  segment  lengths  of  n * 123  sec  and  N = S 
usually  give  good  results. 

The  filters  that  we  call  polarization  and  azimuth  filters  have 
been  applied  before  with  varying  success  at  single  stations  to  shorter- 
period  data  than  we  use  (Archambeau  and  Flinn,  1965;  Flinn,  1965; 

Sax  and  Mims,  1965;  Simons,  1968).  The  major  shortcoming  of  such 
filters,  indeed  of  all  filters,  is  that  they  cannot  improve  the 
signal-to-noi se  ratio  of  the  data  in  any  given  frequency  band.  Their 
strength,  however,  is  that  they  can  objectively  select  and  pass  only 
those  frequency  bands  where  the  signal-to-noise  ratio  is  momentarily 
high.  These  filters,  therefore,  are  .useful  not  for  overcoming  noise, 
which  can  be  done  only  if  more  than  just  single-station  data  are 
available,  but  rather  for  avoiding  noise  wherever  and  whenever  possible. 

Figure  18a  shows  the  application  of  the  Rayleigh  wave  polarization 
filter  to  an  event  from  Venezuela  recorded  at  Ogdensburg.  Comparison 
of  signal-to-noise  ratio  in  the  original  digitized  seismograms  with 
the  signal-to-noise  ratio  in  the  polarization  filtered  seismograms  shows 
an  unreported  event  also  recorded  at  Ogdensburg.  The  signal-to-noise 
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SUHfACE  WAVE  POLARIZATION  WITH  30-60  SEC  BANDPASS  FILTER 
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Fig.  1 8 (a)  An  event  from  Venezuela  recorded  at  Ogdensburg.  (6)  An  unreported  event 
recorded  nt  Ogdensburg.  In  both  (n)  and  (6)  the  vortical  and  north-south  components  arc 
shown.  The  upper  seismogram  of  each  pair  is  the  original  digitized  seismogram;  the  lower 
one  is  the  seismogram  after  polarization  filtering. 
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Fig.  Ha,  The  north-south  component.  The  first  trace  in  the  upper-left  corner  is  the  coda 
of  an  event  from  Iran  (??u  = 6.0,  April  12,  1971,  OT  19h  03m  25.9s)  recorded  at  Ogdensburg. 
T-eneath  it.  is  an  event  from  Honshu  (mt,  = 5.6,  January  4,  1971,  OT  21h  08m  53.4s)  containing 
a Love  wave  recorded  at  Charters  Towers.  This  test,  event  is  busied  with  decreasing  energy 
in  the  coda  in  th?  succeeding  seismograms.  The  seismograms  after  azimuth  filtering  arc 
shown  on  the  right-..  ,.id  side. 
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Fig.  !8b  The  east-west.  component.  The  first,  truer  in  the  upper-left,  corner  is  the  rodn  of 
an  event  from  Iran  (»u>  = f> 0,  April  12,  1971,  OT  lflli  ft'lm  25.9s)  reronled  at,  Ogdensburg.  Be- 
neath it.  is  an  event  from  Honshu  (114  = 5.0,  January  •(,  1971,  01’  21  h 08m  53.-1s)  containing  a 
l.ove  wave  reeordeii  at.  Charters  'Powers.  This  test-  event  is  buried  with  decreasing  energy 
in  the  rodn  in  the  succeeding  seismograms.  The  seismograms  after  azimuth  filtering  arc 
shown  on  the  right-hand  side. 
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improvement  in  the  filtered  records  is  at  least  6 db.  In  the  case  of 
the  unreported  event,  knowledge  of  the  azimuth  was  not  required  by 
the  polarization  filter. 

Figure  19  shows  the  success  of  the  azimuthal  filter  in  extracting 
Love  waves  from  the  coda  of  a large  event.  East-west  and  north- 
south  components  are  shown.  On  the  left-hand  side  of  Figures  19a 
and  19b  are  shown  first  the  coda  of  a surface  wave  and  beneath  it 
a test  event  containing  the  Love  wave  that  we  seek  to  retrieve  after 
filtering.  In  the  succeeding  seismograms,  different  signal-to- 
noise  ratios  are  simulated  by  burying  the  test  event  with  decreasing 
energy  in  the  coda  until  it  is  no  longer  distinguishable.  On  the 
right-hand  side  of  Figures  19a  and  19b  are  the  corresponding  seismo- 
grams after  azimuthal  filtering.  Retrieval  of  the  Love  wave  is 
successful  until  about  the  last  seismogram  from  the  bottom  in  each 
case.  Comparison  of  the  largest  amplitudes  of  the  coda  before  and 
after  azimuthal  filtering  indicates  that  a 14-db  signal  enhancement 
was  achieved. 


2 . Discussion 


Because  both  of  time-varying  filters  have  a built-in  strategy 
of  signal  enhancement,  they  perform  better  than  passive  techniques 
such  as  band-pass  or  matched  filtering.  Our  experience  has  been 
that  the  earthquake  signature  at  periods  longer  than  30  sec  is  not 
really  broad  band  and  hence  makes  an  unsatisfactory  matched  filter. 
By  using  either  the  surface  wave  train  or  the  entire  earthquake  as 
a matciied  filter,  a maximum  s i gna  1 - to-noi  se  gain  of  6 db  is  pre- 
dicted by  the  equivalent  band-width  of  long-period  records,  and 
even  this  is  seldom  realized  in  practice  (Savino  et  al . , 1972b). 

In  some  cases  the  result  scarcely  exceeds  the  s i gnal -to-noi se  ratio 
of  the  original  seismogram.  Band-pass  filtering  often  passively 
rejects  both  signal  and  noise. 

It  is  clear  that  the  success  of  these  filters  depends  on  the 
momentary  content  of  noise.  If  there  is  frequency  separation  be- 
tween signal  and  noise,  the  signal  may  be  extracted  even  if  the  sig- 
nal-to-noise  ratio  is  less  than  1 (Figure  13b).  This  is  an  impor- 
tant criterion  in  trying  to  retrieve  signals  from  a coda  (Figures 
19a  and  19b).  Thus  a long-period  Rayleigh  or  Love  train  that  other- 
wise woul-d  be  above  noise  level  may  be  buried  by  a coda  yet  retriev- 
able at  a single  station. 

If  the  signal  and  noise  share  the  same  frequency  passband,  the 
phase  angle  of  the  signal  is  difficult  to  retrieve  accurately  even 
for  s i gnal -to-noi se  ratios  as  high  as  2;1.  We  would  expect  enhance- 
ment only  of  signals  usually  visible  on  the  original  record.  In 
cases  where  the  signal-to-noise  ratio  is  greater  than  2,  the  filters 
consistently  enhance  the  signal  at  least  6 db.  This  would  still  be 
useful  in  bringing  out  the  character  and  time  of  arrival  of  long- 
period  surface  waves.  However,  these  adaptive  filters  will  cap- 
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italize  on  any  momentarily  advantageous  frequency  bands  when  such 
are  present.  Hence  in  some  cases,  despite  low  signal -to-noi  se 
ratios  and  poor  frequency  separation,  enough  frequencies  may  be 
passed  to  yield  a signal  recognizable  by  its  complex  signature  and 
dispersion.  Unlike  optimum  filtering,  no  a priori  knowledge  of  the 
signal  and  noise  spectra  is  necessary. 

E.  HIGHER  MODES  AS  Afl  AID  TO  SEISMIC  DISCRIMINATION 

A paper  describing  the  application  of  a surface  wave  magnitude  scale 
based  on  the  amplitude  of  higher-mode  Rayleigh  waves  to  the  discrimina- 
tion of  anomalous  earthquakes  from  underground  nuclear  explosions  has 

been  completed  and  accepted  for  publication  In  the  Bulletin  of the  So 

moloqical  Society  of  America  (Forsyth,  1976).  by  using  the  new  higher- 
modescaT  e , I’.’.' , several  anomalous  events  in  central  Asi  . reported  by 
other  investigators  were  shown  to  be  normal  earthquakes.  The  conclusions 
based  on  the  higher  modes  were  confirmed  in  a number  of  cases  by  the 
observation  of  long-period  (40  to  60  s)  fundamental  mode  signals  at  the 
high-gain  long-period  station  at  Chiang  Mai,  Thailand.  A relatively 
high  ratio  of  higher  mode  to  fundamental  mode  energy  can  be  produced 
when  an  unusual  combination  of  source  depth  and  focal  mechanism  leads  to 
poor  excitation  of  the  fundamental  mode  Rayleigh  wave,  in  which  case, 
the  amplitude  of  the  higher  mode  yields  a more  reliable  estimate  of  the 
size  of  the  event  than  the  amplitude  of  the  fundamental  mode.  Several 
anomalous  events  in  central  Asia,  which  are  characterized  by  unusually  . 
low  Ms:mh  ratios  and  hence  could  be  suspected  to  be  nuclear  explosions, 
are  reclassified  as  earthquakes  when  the  new,  hi9her-inode  scale,  is 

employed.  The  focal  depth  and  mechanism  of  these  anomalous  events  ap- 
parently causes  poor  excitation  of  the  fundamental  mode  Rayleigh  wave 
without  significantly  affecting  the  amplitude  of  the  higher  Elim- 

ination of  these  source  related  factors  from  traditional  Ms. mb  discrirn- 
inants  is  an  essential  tool  for  reliable  discrimination  between  earthquakes 
and  nuclear  underground  tests. 


In  compiling  surface  wave  data  for  the  measurement  of  phase  velo- 
cities along  paths  to  the  high-gain  station  at  Chiang  Mai,  great  var- 
iability was  noticed  in  the  amplitude  of  higher  mode  Rayleigh  wave, 
relative  to  the  amplitude  of  the  fundamental  mode.  The  records  from 
the  explosions  and  some  of  the  earthquakes  examined  were  very  s imple  and 
dominated  by  the  fundamental  mode.  Many  other  records  were  quiv.e  complex, 
with  the  amplitude  of  a signal  presumed  to  be  a higher  mode  approaching 
or  exceeding  'the  amplitude  of  the  fundamental.  An  experiment  was  designed 
specifically  to  determine  whether  these  differences  were  source  effects 
or  merely  due  to  changes  in  propagation  path. 

The  effects  of  source  variations  on  the  seismic  signal  can  be  sep- 
arated  from  path  and  receiver  effects  by  examining  sev. eral  closely  spaced 
earthquakes  with  different  source  geometries.  Records  from  three  shallow 
events  (Figure  20  ) are  used  to  demonstrate  the  pronounced  changes  in 
apparent  complexity  of  the  signals  caused  by  relatively  small  variations 
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Trace?  of  digitized  seismograms  observed  at  SUL  on  vertical 
component. 
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Traces  of  digitized  seismograms  observed  at  V00  on  vertical 


component . 


in  source  depth.  Event  2 Is  located  within  about  15  km  of  event  1 and 
event  3 is  approximately  150  km  southwest  of  the  first  pair.  Tracings 
of  the  digitized,  long-period,  vertical  component  from  the  three  events 
as  observed  at  three  World-Hide  Standard  Seismograph  iletwork  (WUSSII) 
stations  are  shown  in  Figures  21,  22,  and  23.  The  three  stations, 

Shill ong  ( S H L ) • Poona  (POO)  and  Shiraz  (SHI),  cover  a 120°  range  of  a z 1 - 
muths  and  are  2,100  to  2,700  kin  from  the  epicenters.  Each  record  has  been 
normalized  to  its  maximum  amplitude  to  facilitate  comparison  of  the  wave- 
forms . 

All  three  records  from  event  1,  top  tracing  in  Figures  21,  22  and 
23,  are  characterized  by  multiple  arrivals  of  short  period  energy  and 
large,  rapid,  amplitude  fluctuations.  In  each  figure,  the  S-wave  is 
near  the  beginning  of  the  record  at  a time  equivalent  to  a group  velo- 
city of  about  4.5  km/s.  At  SHI  (Figure  21),  there  are  arrivals  of 
energy  with  1 js  period  at  3.u,  2.7,  and  2.5  km/s.  There  are  amplitude 
fluctuation  and  apparent  phase  shifts  at  several  points,  such  as  those 
at  about  2.6  and  2.35  km/s.  At  SHL,  (Figure  22)  some  apparent  long- 
period  energy  arrives  at  about  3.0  km/s,  but  is  overidden  by  the  short- 
period  oscillations  which  continue  from  the  S-wave  to  the  end  of  the 
record  shown.  At  POO  (Figure  23),  there  is  a pronounced  arrival  of 
short-period  energy  at  about  3.5  km/s,  followed  by  some  long  period 
signal  (~  30s  period)  and  high-frequency  oscillations  with  large 
amplitude  fluctuations. 

In  contrast  to  the  "messy"  records  shown  at  the  top  of  each  figure, 
some  of  the  Rayleigh  waves  traveling  the  same  paths  from  events  2 and  3 
have  a markedly  different,  "clean"  appearance.  At  SHI  (Figure  21),  the 
record  of  event  2 shows  a very  simple  wave,  normally  dispersed  from  about 
100s  period  at  3.6  km/s  to  about  Ss  at  2.3  km/s.  The  lack  of  apparent 
phase  shifts  and  the  smooth  variation  in  amplitude  do  not  require  the 
path  to  be  homogeneous  between  source  and  receiver,  but  it  does  suggest 
that  it  is  unlikely  that  multipathing  is  primarily  responsible  for  the 
complex  appearance  of  the  record  of  event  1.  Since  the  paths  are  nearly 
i ( II  t i r a 1 , tbCpri!  frOst  - = Stout  fwid  2t,ir_  r.  ta  1 41  f ftrtfitcS  betweeii  the  suut ces 
of  the  two  records. 

In  Figure  24,  contour  plots  of  energy  as  a function  of  frequency  and 
arrival  time  at  SHL,  as  measured  by  moving  window  analysis  (Landisnann 
e t a 1 . , 1 969),  show  that  the  same  dispersion  curve  of  the  fundamental 
mode  Rayleigh  wave  is  found  for  all  three  records  despite  the  distinct 
differences  in  appearance.  The  complexity  of  some  of  the  seismograms 
is  caused  by  .the  superposition  of  relatively  short-period  energy  upon 
the  dispersed,  fundamental  mode  arrival.  Since  much  of  this  energy  ar- 
rives early,  it  cannot  be  scattered  or  reflected  fundamental  mode  energy. 
The  arrival  concentrated  at  3.3  km/s  at  SHL  could  be  interpreted  as  SS 
or  a higher  mode  Rayleigh  wave.  The  nature  of  this  energy  is  much 
clearer  at  PQ0,  particularly  in  the  record  of  event  2,  which  shows  a 
well-dispersed  wave  with  period  about  10s  at  3.9  km/ s to  5s  at  3.4  km/s. 

Trur  disperse  nature  -f  wave  indicates  that  the  signal  i.rhst  be  a 
higher  mode  rather  than  a multiply-reflected  S-wave. 
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There  are  three  primary  cases  in  which  a relatively  high  ratio  of 
hiqhor  mode  to  fundamental  mode  energy  will  be  generated  at  the  source. 

11  The  earthquake  may  be  below  the  depth  at  which  the  fundamental  mode 
is  excited  efficiently,  yet  still  be  sufficiently  shallow  to  excite 
higher  modes  within  the  crustal  wave  guide.  2)  Tor  certain  focal  ro®^n" 
anisms,  fundamental  mode  Rayleigh  waves  will  be  excited  only  weakly  if 
the  earthquake  occurs  at  a depth  corresponding  to  a^zero  crossing  in 
displacement  or  stress  versus  depth  for  the  wave.  Zero  trossings  Tor 
the  higher  modes  are  generally  not  at  the  same  depths,  so  if  the  ampn 
tude  of  the  fundamental  mode  is  reduced  by  this  effect,  the  higher  modes 
may  be  undiminished.  3)  Except  for  sources  very  near  the  surface,  the 
azimuthal  radiation  pattern  for  higher  modes  may  differ  from  the  funda- 
mental modes.  At  azimuths  corresponding  to  a node  in  the  fundamental 
mode  radiation  pattern  there  may  appear  to  be  a relatively  large  amount 
of  higher  mode  energy.  Although  each  of  these  three  mechanisms  for  pi o- 
ducinq  a large  ratio  of  higher  mode  to  fundamental  mode  energy  is 
frequency  dependent,  they  can  affect  a range  of  frequencies  broad  enough 
to  significantly  change  the  appearance  of  the  seismogram.  All  three 
mechanisms  involve  reducing  the  amplitude  of  the  fundamental  mode  rather 
than  increasing  the  absolute  amplitude  of  the  higher  modes.  Thus,  when 
the  higher/fundamental  ratio  is  large  and  the  surface  wave  magnitude 
based  on  fundamental  mode  Rayleigh  wave  amplitudes  is  anomalously  low  _ 
compared  to  the  body  wave  magnitude,  a more  reliable  estimate  of  seismic 
moment  or  surface  wave  magnitude  may  be  obtained  from  the  amplitudes 
of  the  higher  modes. 

As  one  of  the  primary  applications  of  an  improved  magnitude  scale  will 
be  to  aid  in  discrimination  between  earthquakes  and  underground  explosions , 
the  higher  mode  scale  (Mb)  is  designed  to  be  equivalent  to  the  fundamental 
mode  scale  (Ms)  for  a near-surface  focus  event.  The  scale  now  commonly 
accepted  for  Asian  events  (Marshal  1 and  Basham,  1972;  Conference  of  the 
Committee  on  Disarmament,  1972)  is  of  the  form 


Ms  = log  A + O'  (a)  + P(T) 


(1) 


where  A is  the  maximum  amplitude  (nm)  of  the  vertical  component  of  ground 
motion  in  the  Rayleigh  wave  train,  T is  the  period  of  the  wave  of  maximum 
amplitude,  P(T)  is  the  period  correction  which  depends  on  the  propagation 
Dath  and  13'(a)  is  a distance  normalizing  term  that  corrects  for  the  effects 
of  geometric  spreading,  scattering  and  absorption  of  the  propagating  sur- 
face wave.  P(T)  is  assigned  value  0.0  at  the  arbitrary  reference  period 
of  20s. 

To  equalize  M,  and  llh,  it  is  necessary  to  know  the  relative  amplitudes 
of  the  two  modes.  At  the  surface,  the  excitation  function  reduces  to  a 
simple  form  with  only  those  terms  proportional  to  the  displace- 

ment remaining.  The  amplitude  of  the  fundamental  mode  wave  of  20s  period 
should  be  approximately  7.5  times  the  amplitude  of  the  first  higher  mode 
wave  of  10s  period  for  any  very- near-s urf ace  focus  event.  The  formula 
for  Ml?  should  therefore  be: 

Mb  = log  A + B'(a)  + log  T/TQ  + log  7.5 


Fig.  ^6  M,-mb  diagram  summarizing  analysis  of  Tibetan  anomalous  events. 
Arrows  and  diamonds  show  increase  in  Mt  values  when  observations  of  higher- 
mode Rayleigh  waves  are  included,  with  (+)  indicating  original  value.  Arrows 
and  open  circles  show  increase  in  A/s  values  when  observations  of  Love  waves 
are  included.  Dashed  line  contains  all  reported  events  of  the  1968-1969  swarm 
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Magnitudes.  Origin  Timl 

and  Location 

or  Anomalous  Events  and  Explosions* 

Event 

Dale 
(m  a yr) 

Origin  Time 
(li nr.s| 

Liiliuidc 
( Nit  li) 

Longitude 

Depth 

mb 

A 

V; 

05  01  69 

04:00:08.7 

44.0 

77.9 

53 

4.9 

2.9 

3. 

IU 

.17  09  69 

13:41:09.0 

40.1 

70.7 

33 

4.9 

3.6  ,f% 
2.9 

4 

n 

07  01  68 

19: 14: 54  7 

44.0 

79.3 

33 

4.8 

3 

DS 

10  24  71 

08:49:04.6 

28.2 

87.2 

44 

4.8 

2.8 

13 

l-« 

II  24  71 

08:23:24.6 

38.7 

73.3 

33 

4.6 

2.8 

3 

I* 

12  04  71 

08:  38:00.7 

27.9 

87.9 

32 

4.9 

3.2 

4 

Ci 

1 1 30  69 

03:32:57.1 

49.9 

79.0 

0 

61 

4.2 

4 

II 

03  22  71 

04.32:57.7 

49.8 

78.2 

O 

5.8 

4.1 

4 

1 

11  02  72  ' 

01:26:57.5 

49.9 

78.8 

0 

6.2 

3.9 

4. 

* llvpoccmral  data  from  NF.  1C. 

+ All  \/,,  values  have  been  reeomptiled  to  a common  scale  and  may  not  be  identical  lo  original  references. 
} Original  lelercnee:  Landers  ( 19721. 

DOngnial  reference:  Nulllrand  Kmr(1975). 


where  T is  the  period  of  che  maximum  amplitude  A,  B 1 ( A ) is  distance 
correction  from  Marshall  and  liasham  ( 1 972),  and  Tq  is  the  reference 
period  of  10s.  This  formula  should  be  modified  if  used  in  areas  other 
than  central  Asia,  because  O' (a)  and  the  mode  correction  factor, 
log  7.5,  may  vary  from  region  to  region. 


The' ratio  of  surface-to  body-wave  magnitude,  provides  otic 

of  the  principal  means  of  distinguishing  earthquakes  from  underground 
explosions  (e.g.  Evernden,  1969;  Liobermann  and  Pomeroy,  1969).  Explo- 
sions toiid  lv  be  less  ufftetirt  ?tiTi  of  waves*  uitti  M, 

less  than  [m,.  - 1.5].  In  contrast,  IB  for  most  shallow  earthquakes 
is  greater  than  or  equal  to  [m^  - 1.0]  (Conference  Committee  on  Disarm- 
ament, 1972).  However,  there  are  a few  "anomalous"  earthquakes  which 
have  explosion-like  Ms • mg  ratios,  or  fall  between  the  earthquake  and 
explosion  population,  i.e.  [1%  - 1.0]>Me  _>  [m^  - 1.5],  and  thus  are 
statistically  indistinguishable  from  explosions.  Examples  of  anomalous 
earthquakes  in  central  Asia  have  been  presented  by  Landers  (1972),  the 
Conference  of  the  Committee  on  Disarmament  ( 1 972),  Douglas  et  al . ( 1 974) 
and  lluttli  and  Kim  ( 1 975). 


The  long-period,  vertical  component  of  motion  for  two  anomalous 
earthquakes  and  an  underground  explosion  are  shown  in  Figure  25  as  re- 
corded at  the  UUSSU  station  in  Kabul,  Afghanistan  (KBL , Figure  20). 
Locations  and  magnitudes  are  given  in  Table  3.  The  signal  from  event  A, 
an  earthquake  in  Kazakh,  shown  at  the  top  of  the  figure,  happens  to  be 
superimposed  on  the  Rayleigh  wave  train  from  a much  larger  earthquake 
in  the  Tonga  region.  Although  the  high-frequency  Rayleigh  wave  from  the 
Kazakh  event  is  clearly  visible,  it  is  impossible  to  obtain  a reliable 
estimate  of  the  20s  signal  level  of  the  smaller  earthquake.  Event  A 
can  be  seen  at  several  other  stations  in  central  Asia,  but  in  each  case, 
the  surface  waves  from  Tonga  interfere.  All  of  the  anomalous  events  are 
small,  making  it  difficult  to  measure  the  spectrum  over  a period  range 
broad  enough  to  avoid  frequency-dependent  source  effects,  even  when  there 
is  no  interfering  signal.  For  event  A,  Ms  has  to  be  based  on  signals  of 
period  7 to  8s,  yielding  an  estimated  magnitude  of  2.9,  well  within  the 
explosion  population.  However,  as  Landers  (1972)  pointed  out,  the  am- 
plitude of  the  higher  mode  arrival  at  about  3.4  km/s  is  equal  to  or 
greater  than  the  amplitude  of  the  fundamental  mode.  The  revised  esti- 
mate of  surface  wave  magnitude,  til1  = 3.9,  p iaces  the  event  in  the  earth- 
quake population  (Table  3).  Another  similar  event  located  nearby  (C  in 
Table  3)  is  also  removed  from  the  anomalous  category  by  using  Mg.  The 
number  of  observations  on  which  is  based  is  listed  in  the  last  column 
of  Table  3. 


Unlike  events  A and  C,  for  which  no  long-period  surface  waves  could 
be  detected,  event  8 generated  20s  Rayleigh  waves  which  could  easily  be 
measured  at  KBL  (2nd  line  in  Figure  25).  The  M s : ni ^ relation  for  event  B 
falls  between  the  earthquake  and  explosion  populations,  but  again,  the 
amplitude  of  the  higher  mode  restores  it  to  tiie  earthquake  group.  M]  also 
allows  two  of  the  anomalous  earthquakes  listed  by  lluttli  and  Kim  (1975)  to 
be  classified  as  normal  earthquakes  (events  E and  F in  Table  3).  Event  F 
is  unusual  in  that  40-  to  60s  period  Rayleigh  waves  are  clearly  visible 
at  the  UUSSil  station  II D I and  the  high-gain,  long-period  station  at  Clnang 


Mai,  Thailand  (CIIG),  oven  though  very  little  fundamental  mode  Rayleigh 
or  Love  wave  signal  with  period  less  than  30s  can  be  seen  at  any  station. 

Ihe  higher  mode  from  this  event  travels  with  a group  velocity  of  about 
3.B  ktn/s,  in  contrast  to  the  crustal  velocities  observed  for  other  events 
The  mantle-1 ike  velocity  of  the  higher  mode  and  the  long-period  character’ 
of  the  fundamental  node  suggests  that  the  earthquake  occurred  deep  within 
the  imnalaynn  crust,  perhaps  at  a depth  of  70  kn  or  more.  Ils  based  on 
the  oJs  period  waves  is  about  4.0,  in  reasonably  good  agreement  with  IlD. 
Event  D occurred  In  nearly  the  same  location  as  event  F,  but  was  too 
small  to  generate  an  observable  higher  mode.  The  fundamental  node  could, 
however,  be  detected  at  the  CIIG  high-gain  station  where  it  was  recorded  with 
the  same  waveform  at  about  half  the  amplitude  of  event  F.  1 suggest  that 
event  D also  occurred  deep  within  the  Himalayan  crust  and,  by  analogy 
with  event  F,  should  be  assigned  an  !lj.'  of  about  3.8.  Other  examples  of 
higher  modes  from  anomalous  earthquakes  in  the  Tibetan  area  are  presented 
by  Tatham  et  al . (1975)  and  the  effect  of  using  mD  instead  of  Me  is  shown 
in  Figure  26. 

If  the  higher  mode  magnitude  scale  is  correct,  M*}  should  be  equal  to 
Ms  for  very  shallow  earthquakes  or  underground  explosions.  Although  there 
are  as  yet  no  known  explosions  within  an  area  with  a 70  km  thick  crust, 
an  attempt  was  made  to  test  the  validity  of  the  scale  using  signals  from 
the  Soviet  nuclear  test  site  in  eastern  Kazakh  (Figure  20).  One  effect 
of  a thinner  crust  is  to  reduce  the  period  of  the  higher  modes  trapped 
in  ti.e  crustal  wave  guide.  As  can  be  seen  in  the  bottom  trace  in  Figure 
25,  the  maximum  period  of  the  higher  mode  generated  by  the  explosion  is 
only  about  5s.  Mg  for  three  explosions  at  the  site,  events  G,  II  and  I 
in  Table  3,  is  0.05  to  0.15  units  larger  than  Ms.  This  discrepancy  is 
not  large  and  may  be  due  to  the  differences  in  crustal  structure  between 
the  test  site  and  the  locations  of  the  anomalous  earthquakes  for  which 
the  scale  was  designed. 

F.  FOCAL  MECHANISMS  OF  SMALL  MAGNITUDE  EVENTS 

Tatham  and  Savino  (1973)  described  a surface  wave  technique  for  study- 
ing the  focal  mechanisms  of  small  magnitude  events  at  telesclsnic  distances. 
The  method  consists  of  comparing  observed  amplitudes  of  long-period  Ray- 
leigh and  Love  waves  with  radiation  patterns  of  focal  mechanisms  consistent 
with  the  previously  determined  tectonic  regime  of  a certain  region.  It 
Is  important  to  note  that  in  using  this  method  we  do  not  determine  focal 
mechanisms,  per  se,  but  rather  discriminate  between  different  types  of 
mechanisms  th.ought  to  be  prevalent  in  a region. 

In  the  study  of  Tatham  and  Savino  (1973)  the  relative  amplitudes  of 
Rayleigh  and  Love  waves  with  periods  between  10  and  50  seconds  were  deter- 
mined for  some  60  swarm-type  events,  assumed  to  be  shallow  (h<15  km), 
that  occurred  near  the  Delfin  basin  in  the  northern  Gulf  of  California. 

These  surface  waves  were  recorded  at  the  new  high-gain  station  at  Albuquer- 
que, New  Mexico  ( ALQ ) . The  bathymetry  and  seismicity  of  the  northern 
Gulf  of  California  suggest  that  faulting  is  occurring  at  short  spreading 
centers  (dip-slip  faulting)  and  along  transform  faults  (strike-slip  fault- 
ing). The  radiation  patterns  for  Love  and  Rayleigh  waves  corresponding 
to  the  above  faulting  patterns  for  the  northern  Gulf  of  California  are 
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schematically  illustrated  in  Figure  27.  All  of  these  patterns  were 
computed  assuming  double-couple  sources  and  shallow  focal  depths.  The 
position  of  the  nodes  in  the  radiation  patterns  for  strike-slip  faulting 
are  independent  of  focal  depth,  but  there  are  changes  in  the  absolute 
amplitude  of  the  lobes,  especially  for  Rayleigh  waves,  with  variations 
in  focal  depth.  The  dip-slip  radiation  patterns  are  for  a pure  dip- 
slip  mechanism  with  the  fault  plane  15°  from  vertical  and  zero  focal 
depth.  It  is  important  to  note  that  these  dip-slip  patterns  change 
significantly  with  increasing  focal  depth,  even  to  a depth  as  shallow  as 
33  kin.  That  is,  as  depth  increases  for  dip-slip  mechanisms,  the  number 
of  nodes  and  the  number  and  amplitude  of  lobes,  change.  However,  the 
amplitudes  of  the  surface  waves,  especially  Rayleigh  waves,  tend  to 
decrease  as  depth  increases,  an  effect  that  is  especially  pronounced  for 
these  dip-slip  mechanisms.  Further,  the  amplitudes  of  surface  waves 
generated  by  near- verti cal  dip-slip  faulting  tend  to.be  small  at  all 
focal  depths.  The  reader  is  referred  to  the  paper  by  Tatliam  and  Savino 
(1973)  for  more  detailed  arguments  concerning  the  focal  depths  of  the 
swarm  events  included  in  this  study. 

The  near  optimum  azimuthal  position  of  the  high-gain  seismograph 
station  ALQ,  with  respect  to  the  radiation  patterns  shown  in  Figure  27, 
provides  for  the  possibility  of  discriminating  between  the  two  possible 
focal  mechanisms  assumed  for  the  northern  Gulf  of  California.  For 
normal  faulting  on  a spreading  center,  surface  waves  radiated  on  the 
azimuth  of  ALQ  should  be  weak  and  of  approximately  equal  amplitude  dis- 
tribution in  both  Love  and  Rayleigh  waves.  For  strike-slip  mechanisms, 

Love  waves  should  be  of  much  greater  amplitude  than  Rayleigh  waves. 

Figure  28  (a-b)  shows  examples  of  the  relative  excitation  of  Rayleigh 
and  Love  waves  by  two  swarm  events  recorded  on  the  3-component  high- 
gain  seismographs  at  ALQ.  Figures  2Ra  and  28c  show  the  three  orthogonal 
components  of  each  record  segment  normalized  to  a common  magnification  for 
each  event.  It  is  especially  evident  that  the  horizontal  Love-wave  motion 
» dominates  over  the  vertical  motion.  Further,  the  observed  vertical  motion 
does  not  have  the  character  of  fundamental  Rayleigh  waves,  and  may  possibly 
be  due  to  higher  modes.  Figures  28b  and  28d  illustrate  the  result  of  a 
rotation  of  the  horizontal  components  into  equivalent  radial  and  transverse 
traces.  Here  the  Love-wave  motion  of  the  transverse  component  is  clearly 
much  greater  than  the  Rayleigh-type  motion  that  would  be  observed  on  the 
vertical-  and  radial  components.  In  addition,  the  radial  motion  shows  the 
same  lack  of  fundamental  Rayleigh-wave  character  as  that  observed  for  the 
vertical  motion.  These  few  examples  demonstrate  that  Love  waves  are  indeed 
the  dominant  surface  waves  observed  at  ALQ  and  from  this  it  was  concluded 
that  strike-sl.ip  was  probably  the  dominant  faulting  mechanism. 

In  total,  surface  waves  from  GO  events  identified  as  being  part  of 
this  swarm  were  observed  at  ALQ.  In  several  cases,  Love  waves  an  order 
of  magnitude  smaller  than  those  in  Figure  28a  and  b,  an  m^  = 4.2  event, 
were  recorded  demonstrating  the  value  of  the  horizontal  components  and 
, how  they  enhance  the  detection  capabilities  of  a single  seismograph 
station. 
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Pjgure  23  a)  Portions  of  ALQ  sei sinograms  showing  surface  waves  from  two 
swarm  events.  Recordings  start  at  02:3^1,  with  all  components 
normalized  to  the  same  magnification. 

b)  Ihe  same  data  wit.)  horizontal  components  rotated  into 
equivalent  radial  and  transverse  components.  Note  strong  Love- 
wave  motion  on  tx’ansverse  component . 

c)  Same  as  Figure  2a  - for  portion  of  traces  starting  at 
04;f|9. 

d)  Same  as  Figure  2b  - for  portion  of  traces  starting  at 
0M9. 


O'!  <19 


N/S 


E/W 
LOt  LOt 


. I — I 

1IWN 


LOt 


/Avv\/v 


LQt 


04  49 


1MN 


RADIAL 


d. 


V/wvv 


Figure30  a)  VfRA  travel  tine  and  slovness  data.  Squares  indicate 
first  arrivals,  b)  Combined  travel  time  and  slowness  data.  Slope 
lines  plotted  with  slope  of  measured  slowness.  Short  vertical 
lines  indicate  first  arrivals,  c)  First  arrival  data  only,  d) 
Travel  tines  adjusted  so  that  first  arrivals  fall  on  Herrin  et 
al.11  travel  *ime  curve.  {From  Simpson,  Mereu  and  King5) 
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Figure3(  Examples  of  processed  records.  Numbers  over-  array- 
phased-sun  traces  are  observed  slowness  values  (s/deg).  (From 
Simpson,  Mereu  and  King5) 
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G.  HUbSE  DISTORTION  IN  BODY  WAVES 

. .Many  seismic  body  waves  are  associated  with  rays  which  are  not 
minimum  travel  time  paths.  Such  arrivals  contain  pulse  deformation 
due  to  a phase  shift  in  each  frequency  component.  For  sufficiently 
high  frequencies,  the  phase  shift  is  n/2  and  frequency  indeoendent. 
Hence,  the  original  waveform  is  related  to  the  distorted  pulse  via' 
Hilbert  transformation.  The  distorting  effect  of  a f requen cy- i nde- 
pendent  phase  shift  is  successfully  observed  in  seismograms  from  events 
in  several  regions.  The  data  examined  are  long  period  (T>9  seconds) 
and  were  taken  from  the  high-gain  stations  and  the  WWSSI1,  They  in- 
clude deep  earthquakes  (depth>500  km),  in  which  a series  of  well  separ- 
ated S phases  (S,  sS,  SS  and  sSS)  are  available.  These  show  that  the 
waveform  of  SS,  which  .has  been  distorted  in  propagation  through  the 
Earth,  can  be  derived  from  the  waveform  of  sS,  which  is  not  distorted. 
SS  and  sS  form  a Hilbert  transform  pair.  Shallow  events,  in  which 
multiple  S phases  overlap,  also  exhibit  behavior  predicted  by  phase 
distortion.  Rays  critically  reflected  or  refracted  at  a discontinuity 
in  the  Earth  also  suffer  a constant  phase  shift,  which  in  general  can 
have  any  value.  An  important  case  is  SICICS : its  undistorted  wave  form 

resembles  that  of  SKS,  which  has  a minimum  travel  time  path. 

Without  exception,  all  the  distorted  waveforms  bear  little  or  no 
resemblance  to  the  original  waveform.  That  is,  neither  the  first 
arrival  of  energy  nor  the  subsequent  relative  position  of  peaks  and 
troughs  on  a distorted  waveform  appear  at  the  ray  theoretical  times. 
Thus,  r-A  curves  constructed  by  choosing  arrival  times  to  correspond 
to  the  fi  rst  arrival  of  energy  may  be  biased.  Similarly,  doubt  is  cast 
cm  differential  travel  times  chosen  from  first  motions,  or  from  averag- 
ing several  points  on  what  appear  to  be  corresponding  peaks  and  troughs 
0 .two  waveforms.  Some  of  the  rays  most  important  to  seismology,  in' 

nr  d'Stortion  phenomenon  occurs,  include  P and  S (where  d2T/dA>0) 
AB»  SS,  and  SKkS.  Ue  find  that  removal  of  phase  distortion  in 
the  data  is  computationally  straightforward,  and  the  resulting  wave- 
form may  be  exploited  to  full  advantage  in  correctly  picking  arrival 
* i in  6 s a 
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USE  OF  SEISMOGRAPH  ARRAYS  IN  DETERMINING  EARTH  STRUCTURE 
Seismic  arrays  can  be  considered  in  three  modes  of  operation; 


i.cati 

JheJ  r„?Y  »t‘«  «rr«y,  or,  for  a specific  source  location) 

they  can  be  approximated  from  the  calculated  azimuth  and  standard  tables 
of  slowness  as  a function  of  distance.  d lcs 

• Determination  of  seismic  parameters.  Using  signals  from  events 

funrtiWn  1°c®t1°ns»  values  of  azimuth  deviations  and  of  slowness  as  a 
function  of  distance  can  be  determined.  d . 


3)  Event  location.  If  a table  of  slowness  as  a function  of  distance 
is  available,  either  from  2)  above  or  from  standard  tables,  measured  values 
of  slowness  and  azimuth  can  be  used  to  specify  the  azimuth  and  distance  to, 
and  hence  the  location  of,  the  source  of  the  recorded  signal. 

It  is  in  the  second  mode  of  operation  that  arrays  are  able  to 
provide  the  most  useful  information  about  the  internal  structure  of  the 
Ear  tli . 

The  first  distance  derivative  of  the  travel  time  curve  (o'  can  be 
shown  to  equal  the  seismic  ray  parameter  "p"  and  it  is  this  paranuter , as 
a function  of  distance  (see  equation  2)  which  is  required  in  order  to 
derive  a velocity-depth  distribution.  Until  recent  years  dT/da  could  only 
be  obtained  by  smoothing  observed  travel  time  data  and  finding  the  slope 
of  the  smoothed  curve.  With  the  advent  of  modern  seismic  arrays,  it  is 
now  possible  to  measure  dT/da  directly  and  independently  of  absolute  time, 
thus  removing  many  of  the  errors  associated  with  the  earlier  technique. 

The  true  velocity  (v^.)  at  the  maximum  depth  of  • penetrati  on  (dp)  for 
a given  ray  is  related  to  the  apparent  velocity  (va)  at  the  Earth's  sur- 
face (radius  Rjr ) by 

v.  = va  <RE  - dp)  ■ 

t _____  (1) 

It  is  obvious  from  (1)  that  the  slope  (dT/da,  proportional  to  the  inverse 
of  va)  and  hence  the  shape,  of  any  part  of  the  travel  time  curve  is 
determined  by  the  velocity  at  the  maximum  depth  of  penetration  for  the 
rays  arriving  over  that  part  of  the  curve.  Conversely,  if  the  apparent 
velocity  at  the  surface  is  known  as  a function  of  distance,  the  true 
velocity  at  depth  can  be  found  by  determining  the  corresponding  depths 
of  penetration  as  a function  of  distance.  This  is  done  using  the  fol- 
lowing solution  to  the  Abel  integral  equation 

1 n ^E  = I cosh'1  £ dA  (2) 

tkT-tj  > J m 

K o 

where  p = dT/da  or  slowness,  and  the  subscript  1 refers  to  rays  which 

have  reached  a maximum  depth  of  penetration  dp. 


When  the  apparent  velocity  (or  slowness)  at  the  surface  (i.e.  va  or 
dT/da)  is  known  as  a function  of  distance  (a)  the  right-hand-side  of  (2) 
can  be  evaluated,  determining  dp  and,  hence,  from  (1)  the  velocity  at  dp. 
By  evaluating  the  integral  for  various  values  of  a,  different  depths 
are  sampled,  allowing  a velocity-depth  distribution  to  be  determined. 

Figure  29  shows  the  travel  time  curve  and  related  parameters  for  a 
hypothetical  upper  mantle  velocity  model.  Slowness  (dT/da,  in  s/deg), 
reduced  travel  time  (T-a/10  km/s,  in  s)  and  depth  of  penetration  (in  km) 
are  shown  as  a function  of  distance  (a,  in  both  km  and  deg).  Velocity 
(in  km/s)  is  given  as  a function  of  depth  (in  km). 
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This  presentation  stresses  the  importance  of  later  arrivals  in  deter- 
mining velocity  structure.  The  vertical  lines  in  Figure  29  at  the  cross- 
over distances  near  13°  and  24°  show  that  the  first  arrival  portions  of 
the  travel  time  curve  correspond  to  rays  which  reach  their  maximum  depths 
of  penetration  over  a very  limited  range.  Between  depths  of  200  and  700  km, 
first  arrivals  in  this  model  correspond  to  a total  depth  interval  of  only 
100  km.  Depths  near  the  two  discontinuities  at  400  and  650  km  are  always 
related  to  later  arrivals.  Thus,  if  the  details  of  later  arrivals  are 
unknown,  much  of  the  velocity  distribution  is  beyond  the  scope  of  the  data. 


Figure  30  shows  the  complete  travel  time  and  slowness  data  used  in  the 
WRA  upper  mantle  study  of  Simpson,  Mereu  and  King  (1974).  If  only  the 
travel  time  or  slowness  data  in  Fi gure  30a  were  used,  it  would  be  diffi- 
cult to  interpret  any  of  the  later  arrivals.  However,  when  the  travel 
time  and  slowness  data  are  combined,  as  shown  in  Figure  30b  and  d,  later 
arrivals,  especially  between  20°  and  25°  (branch  D in  Figure  31),  clearly 
show  a major  triplication.  When  individual  records  are  considered,  the 
usefulness  of  slowness  information  is  even  more  apparent.  Examples  k to 
o in  Figure  31  show  that  second  arrivals  such  as  that  in  Figure  31k  appear 
as  distinct  phases  only  near  21°-22°.  If  only  travel  times  were  used, 
these  arrivals  could  easily  have  been  misinterpreted  as  belonging  to  the 
later  arrival  extension  of  the  A brand)  (Figure  31,  insert).  As  Figure  30 
indicates,  however,  slowness  determinations  definitely  place  these  arrivals 

on  the  D branch. 


The  basic  information  required  to  derive  velocity  as  a function  of 
depth  within  the  Earth  is  the  slope  of  the  travel  time  curve  (dT/dA,  or 
slowness)  as  a continuous  function  of  distance.  Regions  of  high  velocity 
gradient.,  such  as  the  transitions  in  the  upper  mantle  near  400  km  and  650 
km,  cause  triplications  in  the  travel  time  curve,  so  that  later  arrivals 
must  be  studied  to  obtain  a complete  slowness-distance  relationship.  me 
use  of  seismograph  arrays  to  determine  slowness  not  only  allows  this  para- 
meter to  be  determined  independently  of  absolute  travel  time  but  also  aios 
in  the  picking  of  later  arrivals  and  in  correlating  them  with  particular 
branches  of  the  travel  time  curve. 
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